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, A quantum mechanical wave of a finite size moves like a classical particle and shows 

unique decay probability. Because the wave function evolves according to a Schrodinger 
equation, it preserves the total energy but not the kinetic energy in the intermediate- 
time region of a decay process where those of the parent and daughters overlap. The 
decay rate computed with Fermi's golden rule requires such corrections that varies with 
the distance between the initial and final states, and the energy distribution of the 
daughter is distorted from that of plane waves. The corrections have universal properties 
in relativistically invariant systems and reveal macroscopic quantum phenomena for light 
particles. The implications to precision experiments in beta decays and various radiative 
transitions are presented. 
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1. Introduction: wave zone vs particle zone 

^ : 

lO , A wave length of a particle of a momentum p is given by Planck constant h as h/\p\, 

where h = h/(2w) and is of microscopic size. The momentum eigenstate is a plane wave 
of uniform density and many free waves of a constant kinetic-energy are also uniform in 



o 

£f~) • space and are like free particles. A system of many waves of varying kinetic energy shows 

non-uniform behavior called diffraction. A diffraction pattern has a spatial scale of the wave 
length normally but can become much longer in a system of a space-time symmetry. The 
^ . diffraction of this kind that depends on space-time position in many-body scatterings is 

, studied. 

The diffraction gives corrections to transition probabilities computed by Fermi's golden 
rule. These corrections are connected with calibrations of detectors and might be known 
partly to experimentalists. Even so, it is important and useful to many physicists to clarify 
them. 

In diffraction of light, electron, or others, a potential energy transforms an in-coming wave 
to a sum of waves of different kinetic energies. Now, a many-body interaction transforms a 
many-body state to a sum of the same kinetic energy, and the waves behave like free particles 
and do not show the diffraction at the asymptotic region, t = oo. In non-asymptotic region 
of a finite t, however, the kinetic energy is not constant and takes broad values. So the 
state reveals the diffraction. Since this diffraction is caused by a many-body interaction, the 



TThis author contributed equally to this work 



© The Author(s) 2012. Published by Oxford University Press on behalf of the Physical Society of Japan. 
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 
(http://crcativccommons.Org/liccnscs/by-nc/3.0), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly cited. 




Fig. 1 In two-body scattering, in-coming particles and out-going particles in the particle- 
zone behave like classical particles of constant total kinetic energy but they behave like 
waves of non-constant kinetic energy in the wave-zone. The boundary L$ is the coherence 
length and is normally a microscopic length, but becomes macroscopic in certain situation 
discussed in the present paper. 

pattern has universal properties and appears even in vacuum. The diffraction, furthermore, 
gives peculiar corrections to decay rates that depend on the time interval between those of 
the initial and final states, which we call a finite-size correction. 

Scattering processes are defined with the initial states prepared at t = — oo and the final 
states measured at t = oo, where they do not interact with others and have no interaction 
energy. The initial and final states have the constant kinetic energy and reveal the particle 
nature. Amplitudes and probabilities in the asymptotic region have been studied well [1-5]. 
Now near the scattering center, the states overlap and have finite interaction energy. Thus 
they retain wave natures. We call the former region a particle zone, the latter region a 
wave zone, and a length of boundary a coherence length. Fig. (1) shows them in a two-body 
scattering. In the particle zone, even at a finite t, the states behave like particles. In the 
wave zone, however, the state reveals the wave phenomenon that depends on the position 
and can not be described with only the momentum dependent distribution function [6]. The 
coherence length has been considered microscopic in size of the order of de Broglie wave 
length, which may be true for most cases. Then the phenomena in the wave zone may be 
irrelevant to physics and unimportant. However, there has been no serious investigation on 
this length. We study problems connected ed with the wave zone and find that a new length 
Eh/m 2 , where m and E are a particle's mass and energy, appears for the coherence length 
and becomes much longer than de Broglie wave length in relativistically invariant systems. A 
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space-time dependent phase of a relativistic wave packet (E(p)t — p-x)/h becomes (E(p) — 
p-v)t/h = m 2 t/(hE) of the angular velocity, m 2 /(Eh) at a position moving with the velocity, 
v = p/E(p), as f= vt. The angular velocity becomes small for a light particle or at high 
energy and its inverse gives a new scale of length. The length becomes even macroscopic for 
an extremely light particle such as neutrino. Then the wave zone has a macroscopic size, 
and physical phenomena unique to the quantum mechanical waves occur in the macroscopic 
region. They are natural consequences of Schrodinger equations. The physics in this region 
has not been studied except for neutrino and is the subject of the present work. 

Ordinarily, scattering amplitude is defined in the particle zone and is rigorously formulated 
with wave packets [1, 2], and in practical situations, they are approximated well with plane 
waves. For scattering processes at a finite-time interval, T, in the wave zone, probabilities 
to detect them vary with T and deviates from that of the infinite-time interval. We call 
the deviation a finite-size correction and we study them in various processes involving light 
particles in this paper. 

The finite-size corrections of the scattering amplitude and probability have been considered 
irrelevant to experiments in high-energy regions. Plane waves with a damping factor e _e '*' 
with a positive and infinitesimal e in an interaction Hamiltonian employed often for practical 
calculations are invariant under translations and are extended in space. This method is 
powerful for computing the asymptotic values but does not supply the finite-size corrections. 
Because the amplitude in the wave zone is sensitive to the boundary conditions of the 
initial and final state, that has a dependence on the distance between them. Hence, the 
probability has a finite-size correction which has an origin in the boundary conditions. The 
correction must be included for making a comparison of a theory with an experiment then. 
An amplitude constructed with wave packets implements manifestly the boundary conditions 
and supplies the finite-size correction. 

Previous studies of decay processes at finite-time intervals in the particle zone using an 
interaction Hamiltonian of the damping factor e _e '*' [7-10] showed that time dependences 
of decay law of unstable particles are modified from simple exponential behaviors due to 
higher-order effects. These analysis and others of computing the decay rates are applicable to 
kinematical regions where the wave functions of the parent and daughters do not overlap. As 
was pointed out in Ref. [8] correctly, the standard method can not be applied in kinematical 
regions where they overlap. The states have the wave natures, and the decay rate and other 
physical quantities in this region have been thought neither meaningful nor computable since 
then. This is the region, in fact, where the probability to detect the decay product has a large 
finite-size correction. It is a main subject of the present work to develop a theory of S-matrix 
that satisfies the boundary conditions of the measuring processes and to find formulas of 
the physical quantities in this region. One of our results for decay rate T(T) at the large 
distance L = cT, (T < r) is 

T(T, er) = r + N 7f^2 F2 {-m 2 ), m 2 = m 2 parent - m 2 daughter , (1) 

where To is the asymptotic value, r, a, E, m, rridaughter, and m pa rent are the mean-life time, 
wave packet size, energy, and mass of the detected particle, and the mass of daughter and 
parent, respectively. N is a numerical constant and F is the form factor. The second term 
in the right-hand side of Eq. (1) is inversely proportional to T and vanishes at T — > oo. So 
this is the finite-size correction. From its form, the correction becomes significant in small 



3/57 



m, large a and E, and become macroscopic in their sizes for light particles such as photon 
or neutrino. This shows 

hm { hm r(T>)) =r , (2) 

lim ( lim r(T»l = oo. (3) 

T->oo lcr->oo J 

In Eq. (2), the rate becomes the asymptotic value, whereas in Eq. (3), the rate diverges. The 
energy distribution reveals also unusual properties even at T — > oo, if particles of large and 
small sizes are involved in one process. They should appear in various situations such as an 
interface of two phases, and interesting physics is expected. Implications on particle decays 
are studied. 

The transition probability P composed of many processes in the particle zone are factorized 
to that of each microscopic process, Pi, as 



p = U p > ( 4 ) 



Now, the probability for transition processes in the wave zone is not factorized due to the 
finite-size corrections, but the whole process is described by the product of wave functions 
of each microscopic process 

i i 

Because the probability of the whole processes is not factorized, Markov nature of the mul- 
tiple processes is lost. The non-Markov nature is related with a EPR correlation [11] and 
may give various implications. 

The decay rates are studied in the present paper and the scattering cross sections are 
studied in a subsequent paper. This paper is organized in the following manner. In Section 
2, a wave function and S-matrix at a finite-time interval are shown to be different from 
those of the infinite-time interval. Particles described by wave packets and their interactions 
caused by a local Hamiltonian are summarized in Section 3. Two-body decays are studied 
in Section 4, and radiative decays of atoms, nuclears, and particles are studied in Section 5. 
In Section 6 we study decay processes and thermodynamics of quantum particles. Summary 
is given in Section 7. 

2. A finite-time interval effect 

In a physical system described by a Hamiltonian H composed of a free term Hq and an 
interaction term Hi n t 

H = Hq + Hi n t, (6) 

a wave function \^(t)) follows the Schrodinger equation 

ih^(t)) = (H + H int )\^(t)). (7) 

In field theory, the free part Hq is a bi-linear form of fields and the interaction part Hi n t is 
a higher polynomial of fields. Hi n t causes a change of a particle number such as a decay of 
a pion into a charged lepton and a neutrino. 
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2.1. Finite-size correction to Fermi's Golden rule 

A transition rate from an eigenstate of Hq, \a) of energy E a , to another |/3) of energy Ep, 
in a wave zone at a finite-time interval T, seems to be computed with the amplitude / and 
probability P [12] in the form, 

/= C dt(p\H int (t)\a) = [ T dte'^-^F^, (8) 
J o J o 

F ajfj = (p\H 1 (0)\a), 

P = \F a ^\ 2 D(Ep-E Q ;T), (9) 

nfF F . T) _ 4sin 2 [(£ /3 -£ a )T/2] 
D(Ep - E a , T) ^- _ , 

where -F a)J g is the matrix element. In a particle decay, the final state constitutes two or 
more particles of a continuous energy spectrum and the oscillating function D{Ep — E a ;T) 
approximately agrees with Dirac's delta function at the infinite T [13-15], 

D(Ep - E a ; T) = 27rT5(Ep - E a ). (10) 

Because the integral of a function F(Ep) with the weight D{Ep — E a ;T) over the energy 
E/s is computed with a variable x = (Ep — E a )T as 

P = [ Ea+AE dEpD(Ep - E a ; T)F(Ep) 
JE a -A E 

= T / a ' T dl fW 2 F(l/T) , 



J ~A E T \ J 
F{Ep) = \F a ,p\ 2 . 

The symmetric region of the integration was chosen in Eq. (11). At a large T, F(x/T) is 
replaced with F(0), and Eq. (11) becomes 

P = TF(0) / AET dx ( Sm(x/2) ^ 2 = 2^TF(0). (12) 
J~A E T \ X J 

Thus the transition probability integrated over final states is given by 

P = 2ttT J df35(E a - Ep)\F a ,p\ 2 , (13) 

and the rate P/T is constant. This is Fermi's golden rule. 
Now at a finite T, expanding F(x/T) in a power series of x/T 

F(x/T) = Y J C l (^)\ (14) 
l 

we have Eq. (11) in the form, 



p = yjqTi-</^4^p. (i 5 ) 



The integrals over x is easily evaluated. In a small \x\ region, the integrand vanishes for I > 1 
and is constant for I = 0. In a large x region, the integrand behaves as ^x l ~ 2 . So the above 
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integrals become 



2vrTC + C l T 



l>i 



r-A E T 
I dx 

-A E T 



.r 



1-2 



2 



I - 1 



i-i 



2ttTC + °i T 



l>2 



27rTCo < 1 + J] 



l>2 



Q Ae 1 - 1 
Cq 2vrT(/ - 1) 



(16) 



1/T correction is in the second term of the right-hand side, which is finite if A^; is finite. 
The 1/T correction depends on Ag and the distribution of eigenvalue and converges if Ae is 
finite. Appendix A and B study 1/T of various distributions. The value at T — > oo is defined 
uniquely then. 

In relativistically invariant systems, A^ = oo and the correction for I > 2 in Eq. (16) 
diverges. The infinite correction emerges due to a large overlap of wave functions in the 
situation where the ordinary scattering theory can not be applied [8]. The probability at 
a finite time measured with an apparatus does not diverge. Hence the amplitude defined 
according to boundary conditions of measurement process should give the finite value. The 
boundary condition at the finite T is different from that at the infinite T, hence the ampli- 
tude that satisfies the boundary condition at the finite T is different from those of the 
T = oo. As is seen explicitly later, the function decreases rapidly with x/T in such manner 
as e~ a ( x l^ 2 in the right-hand side of Eq. (14), where a is a size of wave functions determined 
from the boundary condition, and the coefficients converges. Since the amplitude at large 
x/T is determined by the boundary condition, the 1/T correction becomes a finite value 
that depends on the boundary condition. Nevertheless they follow a universal relation. It is 
a subject to find the universal properties of the finite-size corrections. 

The states \/3) of satisfying Er = E a contribute to the decay rate, Eq. (13) and the states 
\(3) of Ep 7^ E a contribute to the finite-size correction. Since Ep is continuous, those states 
of Ea ~ E a are sensitive to boundary conditions and so is the finite-size correction. For the 
computation of the probabilities of processes measured in experiments, the wave functions 
for the outgoing waves and in-coming waves should be localized around their centers, as 
were emphasized in textbooks of quantum field theory [16]. The wave packets satisfy this 
property and are used. They can be replaced with the plane waves in the particle zone, but in 
wave zone they are necessary to compute the finite-size corrections. We study the finite-size 
corrections to transition probabilities with wave packets. These amplitudes are described by 
scattering matrix S[T] defined at a finite-time interval T. S[T] is different from the ordinary 
S-matrix, 5[oo], and has unique properties. The coefficients in 1/T expansions become finite 
with a use of S[T]. 

2.2. Wave function at a finite time 

S[T] satisfies the boundary condition at T, and is determined by wave function which starts 
from an initial state at t = and ends at a final state at t = T. Kinetic energy is not a good 
quantum number in the wave function at a finite T. A time-dependent solution of Eq. (7) in 
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the first order of H%nt that satisfies a initial condition 



|tf(0)> = |^ (0) ), H \^)=E \^), (17) 



is 



|*(t)) = e-^ ||V (0) ) + J dpd(u,t)\f3)(P\Hint\^ 0) ) \ , (18) 
u = Ep-E , H \p)=E f) \p), 



where 



d(w,t) = = -2i — K —L-L e —<* m ( 19 ) 

At i — > oo, d(u>,t) becomes 

d(w,t) =2mS(u), (20) 

and the wave function 

Mt)oo) = {|^ {0) > + 27Ti\f3)(f3\K m t\^ 0) )\E,=E } ■ (21) 

has the kinetic energy, Ep ^ Eq. At a finite i, on the other hand, Eq. (20) is not fulfilled and 
the wave function is a superposition of wide spectrum of the kinetic energy, Eg. An average 
of d(u,t) over a finite-time interval St of satisfying ui5t 3> 1 is 

d(u>,t) = , (22) 

and the average of the wave function over the finite interval is 

|*(i)avera g e> = ^ ||^) - 2ttz J d^\0) (/3|F int |^°>) \ Ef) =E } ■ (23) 

In both cases, the state vectors (<)«>) an d ^(i) average) have the frequency Eq/H and the 
total energy Eq, 



H\V(t) oo ) = E \y(t) oo ), (24) 

ff|#(i)average) = Eq |* (*) a verage) • (25) 

Thus the wave function at a finite time £ is a sum of those of broad energy spectrum of Hq 
, whereas that is composed of a discrete spectrum Ep = Eq at t = oo. The conservation law 
of energy defined with H is reduced to the conservation law of the kinetic energy defined by 
Hq only at t = oo. 
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2.3. Scattering operator at a finite-time interval 

Physical quantities are observed through scattering or decay processes and are computed 
with S[T], which is defined from unitary operators 

U(t) = e~ iHt , U = e~ iHot . (26) 

M0ller operators are defined in the form 

Ot(T)= lim UitfU^it), (27) 

and satisfies 

e iHet n T (T) = n T (T ± e t )e iH ° et . (28) 
Scattering operator at the finite time T is the product 

S(T) = nL(T)n + (T), (29) 

and satisfies 

S(T)H = H S(T) + i|An_(T)J n+(T) - inl(T) A n+ (T), (30) 
and the commutation relation 

[5-(T),fl-o] = i|^n_(T)| n+(T)-inL(T)^n+(T). (si) 

Thus S(T) does not commute with Hq, and the conservation law of kinetic energy is violated 
at a finite T. 

From Eq. (31), the matrix element of S(T) between the eigenstates of Hq has energy 
conserving and non-conserving terms, 

(fi\S(T)\a) = 5 e (E a - Ep)f{T) a> p + 6f. (32) 

where the second term, 5f, vanishes at the energy Ep = E a . Since the energy Ep of the first 
and second terms are different, the total transition probability is a sum of each probability. 
The first term gives a normal constant probability that is computable also by the ordinary 
S-matrix of plane waves, whereas the second term gives a T-dependent correction that is 
not computable by the ordinary S-matrix. In ordinary situations, the energy non-conserving 
terms are negligible but they are important in situations of the present work. 

A magnitude of 5f and a probability derived from 5f depends on a dynamics of the system. 
When E a and Ep are approximate energies of the states \a) and |/3), we have 

(E a - Ep)<J3\S(T)\a) = W\0(T)\a), (33) 

o(t) = J A _(T)| n + (T) - inL^n+(T). 



d .. ...V.. _ ._ t ^ d 



Hence 



Sf= mim, ,34) 

&a — 
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and the transition probability to the energy non-conserving states is given in the form, 

where the equality is satisfied at T — > oo. States at ultraviolet energy regions couple in a 
universal manner with the operator 0(T) and contribute to the probability at the finite- 
time interval. Since states of unlimited-momentum couple in Lorentz invariant manner, they 
give a universal correction to Eq. (32). The finite-size correction appears even in the lowest 
order of perturbative expansions and is useful for probing the physical system at the large 
momentum region. 

Boundary conditions necessary to determine a solution of a wave equation uniquely in 
scattering or decay processes are asymptotic boundary conditions [1]. For a scattering from 
an initial state \a) to a final state |/3) of a scalar field expressed by 4>(x), the states \a) at t = 
—T/2 are constructed with free waves 4>i n (x) and the states at t = T/2 are constructed 
with free waves (fc ut(x) an d satisfy asymptotic boundary conditions, 

Urn AaWitW) = H^Ll£> ; (36) 

t->— T/2 

lim (a\<f/(t)\P) = <o#Ll/3>- (37) 
The expansion coefficient $ (t) is defined by 

(f)f(t)=i I d 3 xf*(x,t)d o 0(x,t). (38) 



(f){ n and <$mit are defined in the same way. C-number functions f(x,t) are normalized and 
satisfy a free wave equation. The normalized functions decrease fast in space and form a 
complete set with those functions translated in space. Hence they have central values of 
position and momentum and the state vector is specified by both variables as \p, X). Thus 
matrix elements of S[T] are defined as (pi, Xi\S[T]\pj, Xj) and depend on the position and 
momentum. The finite-size correction are computed with the position dependence of the 
probability. For normalized functions to form the complete set, those of different center 
positions are required [17]. Those of the initial state represent the beam and those of the 
final state represent a detected particle. They are determined by the experimental apparatus 
and those of the initial and final states are different normally. Being non-normalizable, plane 
waves are not suitable for these functions if the damping factor e _<E '*' is not included. Instead, 
wave packets are normalizable and are suitable. <j)in(x) and 4> ut( x ) satisfy the free wave 
equation and the states \a) \f3) are defined with wave packets. The wave packets which have 
finite-spatial sizes and decrease fast at large \x — xq\ ensure the asymptotic conditions at 
a finite T, where xq is the center's position. Hence S[T] is described by wave packets and 
the finite-size corrections are studied with S[T]. We present several examples where the 
finite-size corrections are non-negligible and give interesting observable effects. 

3. Quantum particles described by wave packets 

Waves of finite sizes expressed by wave packets used for formulating S[T] exist in various 
areas. Wave functions at the particle zone lose wave nature quickly and time evolution of 
objects turns to be described by classical equation of motion. Thus a classical mechanical 
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description is smoothly obtained starting from the quantum mechanical description and 
physics in this region is understood well by classical physics, as is explained in most textbooks 
of quantum mechanics. Now in the wave zone, the phase of the wave function remains 
and gives physical effects that are different from classical physics. They could appear in 
macroscopic space-time regions. Then their universal properties are common in any wave 
functions, and can be studied with Gaussian wave packets. Their sizes are determined from 
physical processes of the particles. 

The physics of quantum particles has been neither completely explored nor understood 
and is becoming relevant to recent advanced science and technology, especially for precision 
experiments of light particles. Various phenomena of neutrinos and photons caused by these 
unique phases are studied hereafter. Neutrino interacts extremely weakly with matter and is 
not disturbed by environment, hence its phase is not washed out, consequently the neutrino 
retains the wave nature even in macroscopic area and reveals large finite-size corrections 
[18, 19]. The finite-size correction is observed as a diffraction pattern of the neutrino produced 
in a decay of the pion and in other processes that the neutrino gives rise to. 

Photon is massless in vacuum and behaves approximately like a particle of small mass 
in high-energy region in dilute matter. Normally, the quantum mechanical phase of single 
low-energy photon is washed out and a large number of these photons behave like a classical 
electromagnetic wave in macroscopic areas. In various exceptional situations, its phase is not 
washed out and photons reveal unusual properties and interact with microscopic objects as 
single quanta. The photon then is expressed by a wave packet and the probability to detect 
it at a finite distance shows the diffraction behavior of a single quantum. They lead also the 
photons to have unusual thermodynamic property. 

Waves of small sizes move like classical particles [20-30] and exhibit wave-like behaviors 
such as anomalous finite-size corrections in scattering cross sections or decay rates, and are 
called quantum particles in the present paper. Quantum particles of relativistic waves have 
universal properties. 

3.1. Symmetric wave packets 

The Gaussian wave packet of a relativistic particle of mass m and the central momentum 
Po, position Xq, and time To is expressed in the momentum representation by 

(t,p\p , X , T ) = jV CT 3/2 e -i£(p)(t-T )-^.X -f {V~P )\ (39) 

where a is a spatial size of the wave packet, N is normalization factor and the energy is 
given by a relativistic form, E{p) = \/p 2 + m? . This is a superposition of eigenstates of the 
energy and momentum of the widths an d ^= respectively, and is a simple Gaussian 

form of p at t = Tq, and retains its shape afterward. The completeness of wave packets of 
the continuous center position and momentum and other important properties were given 
in Ref. [17]. Some of them are summarized in the following for a completeness of the present 
paper. They satisfy 

/ ^7^3^^ T >(^^ T I = !. ( 4 °) 
and the wave function in coordinate representation is 

w(p ,x) = (t,x\p ,X ,T ) = dk(x\k}{t,k\po,Xo,T ), (41) 



10/57 



and becomes also a Gaussian form in x around a new center 

u;(pb,2;) = Ar e -^( £ -^-^(*- T «)) 2 e- iS ^ )('- To ) +i P°-^-^ ), (42) 

in a small |i — To| region. Thus the wave function keeps the shape and moves with a velocity 
vq and the modulus is invariant under 

t-tt + dt, x -> sc + w 5t. (43) 

Since the position of the wave packet moves uniformly with the velocity vq and has the 
extension a, the wave function becomes finite only inside a narrow strip of this width. Hence 
the quantum state expressed by this wave packet behaves like a particle of the extension a. 
At a large \t — To|, the function expands. 

The wave functions Eq. (42) decrease rapidly with \x — Xq — vo(t — Tq)| and vanish at 
\x — Xq — vo(t — To) I — > oo. Hence they satisfy the asymptotic boundary conditions and are 
appropriate to use as basis, f(x, t), of Eq. (38). The transition process of the particle prepared 
at the initial time Tj and of observing the final states at a final time Tf of a finite T = 
Tf — Tj, is studied with the S-matrix at the finite-time interval S[T] thus defined. Because 
the S-matrix of plane waves defined at T = oo, <S[oo], satisfy the boundary condition at 
t = ±oo, that is different from S[T] defined at t = ±T/2. S[T] defined by the wave packets 
Eq. (42), and the amplitudes and probabilities obtained from them are not equivalent to 
those obtained from S[oo] generally in the wave zone. Then, the computations should be 
made with S[T]. Conversely if they are equivalent, the computations can be made with either 
methods. 

The kinetic energy is strictly conserved in both classical collisions of particles under a force 
of finite range and quantum collisions described by S[oo] of the stationary states of the free 
Hamiltonian, whereas the conservation law of the kinetic energy is slightly modified in a 
collision of the finite-time interval T described by S[T] from the algebra Eq. (31). The total 
energy is conserved, but is different from the kinetic energy in the space-time region where 
the interaction Hamiltonian has a finite expectation value. Hence the kinetic energy is not 
conserved in this region. The non-conservation of the kinetic energy is a unique property of 
quantum particles described by S[T] and causes unusual behaviors to the collision or decay 
probabilities. 

The quantum states of finite-spatial extensions are expressed by superpositions of plane 
waves of different momenta and energies, and their scatterings are those of the non-stationary 
states. These non-stationary wave packets are specified by the values of position, momentum, 
and complex phase at the center. Even though its spatial size is so small that it behaves like 
a point particle, a wave nature represented by the phase remains. The phase that depends 
on dynamical variables gives physical effects that are characteristic of the quantum particles. 

a in Gaussian wave packet determines the spatial size of quantum particle, and depends 
on the situations. Because the probability to detect this particle is unity inside the wave 
packet, this size is a classical size of a quantum particle. So, a for out-going state is the size 
of the unit of detecting system that gives a signal, and is a size of nucleus used in detector 
for the neutrino. For a high-energy photon, the signal is taken from its e + e~ creation around 
an electric field of the nucleus used in detector, hence a is about the size of nucleus, a for 
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in-state is also the size of wave function that expresses this particle. This size is infinite 
for the idealistic particle in vacuum, but is finite in matter due to effects of environment. 
When a particle expressed by a certain wave function interacts with others and the both 
make a transition to other states, this particle is expressed by one wave function in a finite 
time-interval between these reactions. Hence that is determined by a mean free time of this 
particle. Thus a is determined by the mean free path for in-coming waves, a for pion, kaon, 
muon, proton, photon, and electron in the initial states are estimated from their mean free 
paths in the matter of experiments. Actually most of them are macroscopic sizes in high- 
energy regions. An electron easily loses energy by electromagnetic showers and is exceptional. 
In low-energy regions, an electron, negative muon, and negative pion form bound states 
of microscopic sizes with nucleus in matter, and a are microscopic sizes. Positive charged 
particles such as positive muon and positive pion do not form the bound states with nucleus 
and may have larger a. 

Thus a of a nuclear size, atomic size , or larger size appear depending on the situation. In 
a scattering or decay of waves with different sizes, the wave functions overlap in finite and 
asymmetric region. Consequently conservation laws derived from space-time symmetry are 
modified. 



3.2. Local interaction 

Characteristic features of quantum particles are connected with the phase factor of wave 
functions and appear in the lowest order of interactions of scalar fields. Hence we study 
scattering of particles caused by the local interaction 

j=N 

J?int=gil<p j (x), (44) 
3=1 

in the lowest order of g first. Effects of spin and internal structure will be included later. 
Interactions of Ai in-coming and A2 out-going particles expressed by the wave packets 
parameterized by (pi, Aj, Tj; <7j) at a space-time position (t,x) is given in the form [17], 

j=N ATi N 2 

01 n Vj( x )\ l ) = Y[w* k (x,p k ;X k ,T k ,a k ) x JJ w t (x,pi; X u T h 07) 

j=l k=l 1=1 

= Ar te -^7( a? -^W) 2 -^7(*-*o) 2 e fi+i0 ! (45) 

N t = HiqN lt 

k,l 

where a s and at in the exponent display the extents in x and t and are expressed in the 
form: 

± = E ±,± = E S_5?, (46) 

a s a,- a t aj a s 

j j 

3 
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Here, xo(t) is the center in x and moves with vo, 

x (i) =v o t + x {0), (48) 

^o(O) = a J J2 ^ - *£(±)# 



J J 

X 3 = X 3 ~ VjTj. 

In above equations and hereafter, (+) and (— ) are for in-coming and out-going states, 
respectively. The real part of the exponent of Eq. (45), R, determines the magnitude 
and is composed of the position-dependent and momentum-dependent terms. The former, 
Rtrajectory, and the latter, Rmomentum, are expressed by, 

R = Rtrajectory + Rmomentum, 

JW-, = - £ § + (e + f £ ^S^) • < 5 <» 



2cr 7 ' \ ^ — ' 2<j-j } \ 2<jj 

3 \ 3 



P - 

-^momentum — ^ 



2 



]T (±) (E®) - } - ^ ( I ■ ( 51 ) 



From Rtrajectory, the particles follow classical orbits and from Rmomentum, Eq. (51), they 
follow the approximate energy-momentum conservation. Because the interacting system is 
invariant under a translation of the coordinate system, Rtrajectory is invariant under the 
translation 

Xi -> X % + £ Ti -> Ti + 5, (52) 

where (6,d) is a constant four vector. From Rmomentum, the momentum is approximately 
conserved with the uncertainty 1/^/ctJ and the energy of the system moving with vq is 
approximately conserved with the uncertainty 1 / y/at. Since a massless particle has the max- 
imum speed, the moving frame has a large velocity and the effect becomes significant for the 
massless or extremely light particle. The product Eq. (45) depends also on the phase factor 

= 00 + 01) (53) 
<j> = J2^)(PrX j -E(f j )T j ), 



-2* t E 



(v - Vj)-Xj 



2a j 



where 0o agrees with that of a plane wave. 
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When the values of a s and at are finite, the product Eq. (45) becomes finite in a small 
region of (t, x) and decreases steeply away from this region. Hence the integration over (t, x) 
becomes 

j=N 

J d 4 x{j\ H <fj(x)\l) = N t (2a s 7r)l(2a t )h R+i t (54) 
j'=i 

and converges fast. The integral over < t < T becomes 0(e 2<Tt ). Thus the finite-size cor- 
rection to the probability is 0(e 2at ) with a microscopic a t , and is negligible at a macroscopic 
T. 

3.3. Pseudo-Doppler effect 

First effect caused by the modified conservation law of kinetic energy is the distortion of the 
energy distribution, which appears in the amplitude at a finite and infinite T. 
The energy-momentum conservation in an invariant system under the translation 

Xa ->■ x a + dp, (55) 

where d a is a constant four vector, is derived from the integration for the plane waves 

' dW(**- fc '>* = (27r)V 4 )(^ - k f ), (56) 

where fe, and kf are the four-dimensional momenta of the initial and final states. In the 
amplitude of the wave packets, the wave functions overlap in a finite space-time area and 
the amplitude is not invariant under Eq. (55), but is approximately invariant for a small 
at under the transformation Eq. (43). Hence the energy in the moving frame Eq. (51), is 
approximately conserved. In a system of at = oo, the invariance is rigorous, 
^momentum is rewritten as 

^momentum = ~^(SE) 2 - y(^) 2 , (57) 

5E= E (Ef(p k )-v .pf), 

initial, j final, k 



For small a s and large at, \Sp\ becomes large but \5E\ becomes small, and the modified 
conservation law, 5E = 0, 

E ^Vi) - E Ef (pk) = E • pi - E • ( 58 ) 

final initial final initial 

is fulfilled. The momentum spreading is large and the conservation law for the events of 
5p = or vq = becomes the form: 

E E\p 3 ) - E E f (Pk) = 0. (59) 

final initial 

For the events of 5p ^ and i/o ^ 0, the law becomes 

E V&<Pi) - E 7kE f (p k ) = 0, 7i = ^E^ ^^ (60) 

final initial ^ 
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where jj is the rate of the energies in the moving and rest systems. 5E is also written in 
high-energy region in the following form: 



From the momenta and energies of particles in the final state, f)j can be computed from 
Eq. (61), and E(pj) is calculated. Then Eq. (58) can be verified. The total momenta are 
distributed with the width given by 1/ 1 ^Ja~l but the sum of total energies at pj vanishes 
at each event. Even though the detector is at rest and a real Doppler effect is irrelevant, 
the kinetic energy of the moving frame, instead of that in the rest system, is conserved. 
Consequently the kinetic energy of the final state shifts in magnitude in events of large \Sp\. 
In Doppler effect, the energy shifts in all events, so the shifts due to the wave packet is 
different and called pseudo-Doppler effect. 

For small a s and at, \5p \ and \5E(pj)\ become large. For plane waves, cr, = oo, the velocity 
vq vanishes and the modified conservation law becomes the standard one. Thus the energy- 
conservation for the wave packets is different from both of the classical mechanics and of the 
plane waves of the quantum mechanics. 

The modified law of energy conservation results from S[oo] that satisfies the commutation 
relation [£[oo], Ho] = due to the fact that the wave packets are superpositions of states of 
continuous eigenvalues of H$. The quantum particle of the momentum and kinetic energy, p 
and E(p), and size a makes a reaction as a particle of the energy ^yE(p) and the modified 
conservation law, Eq. (60), is fulfilled. Here ji is regarded as a ratio of the time intervals in 
the moving and rest frames and Eq. (58), is understood as that for the average values taken 
over the time intervals 



Thus the conservation law of energy is modified to that for the average values. Because the 
energy is conjugate to the time, the equality of average values taken over the time intervals 
is reasonable. From Eq. (58), the effective action 



of the initial state coincides with that of the finial state in the present reaction. 
3.4- Finite-size correction 

Second effect caused by the modified conservation law is the large finite-size correction. If 
at is finite of a microscopic size, the integration over t converges and the amplitude and 
probability decrease rapidly due to Rtrajectory In a marginal case of at = oo, modulus of 
Eq. (45) does not decrease with t but wave packets overlap in the infinite-time interval. This 
happens in various situations. If all the particles except a particle 1 are plane waves, 




(61) 




(62) 




(63) 



a x ± oo, 

Oi = oo, i ^ 1 



(64) 
(65) 
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the frequency and the real and imaginary parts of the amplitude 

II N j II N^-^-^^-^-^e 11 ^, (66) 



j I 

a s = ai, a t = oo, 



are 

UJ 



Ei-J2 E J-vi- Sp, (67) 

initial final 



R — ^momentum ^trajectory j 
Rmomentum — ^"(^P) ! ^trajectory — 0, 

cj> = E(p 1 )T 1 -p 1 -X 1 -5p-X 1 . 

The modulus of Eq. (66) decreases fast with \x — xo(t)\ and the total momentum is approxi- 
mately conserved, whereas that is constant in t. The phase factor has a similar form to that 
of the plane wave Eq. (42) but the angular velocity uj is not identical, uj in Eq. (67) is an 
energy in a moving frame of the velocity v\ of showing the pseudo-Doppler effect, and is 



u) = -.Ei (pi) +vi-pi + uj = -(yPi+mf - |pi|) +w 

= -^,+"0, (68) 

where rrt\ is the mass of 1 and ujq is independent of pi at high-energy region. Hence uj depends 
on the momentum of particle 1 differently from the plane wave, and there are more states 
satisfying uj ~ than those of the simple plane wave. The amplitude Eq. (66) at a large-time 
interval is determined by a state of satisfying ui = and also the states of uj ~ 0. From 
Eq. (68), uj is degenerate at |pi| — > oo, and infinite number of states give the contribution. 
It will be shown that the rate derived from this at a finite T, T(T), has the large finite-size 
correction and is described in the form, 

r(T) = r + r 1 (T), (69) 

ri(T) = Ci/T, 

where C\ is a constant and Tq is the asymptotic term. 

at becomes infinite when the right-hand side of Eq. (46) vanishes. This condition is fulfilled 
in particular momenta of the initial and final states. The probability has the finite-size 
correction in this kinematical region, then. 

3.5. Asymmetric wave packet 

In some situations, the wave packet is asymmetric in ki and kj- which are parallel and 
perpendicular to the central momentum, or in k and |fc|. The small energy uncertainty, 
5E <C \k\ also often appears. For an asymmetric wave packet or a wave packet with different 
spreadings in the momentum and energy, we have 

<i,p1p-o,*o,ToU = ATafe-^^)^)-^^-^^-^) 2 -^ (ftO", (70) 

(t,p\p ,X ,T ) E = iVafe-^^^-To)-^-!^) 2 -^^^)-^)^ (n) 

where ctl, &t, and cte are a size in the parallel and perpendicular directions to the center 
of momentum, and that in the energy. The functions in coordinate representation become 
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Gaussian forms in x and t. 



{t,x\po,Xo,T ) asy = J dp{x\p){t,p\p ,X ,T ) asy 



= iVcrf e - iE (^)( t -^o)-ip a -X -^{x L -X L -v{t-To)f-^{x T ) 2 ^ ^ 

(t,x\po,X ,T ) E = J dp (x\p)(t,p\p ,X , T ) E . (73) 

In ot or oe ~ c, the energy spreading is about the same as that of the momentum, and 
the probability for a finite |Ap| around the central momentum po has the pseudo-Doppler 
and finite-size effects. In ai S> &t or oe 3> a, on the other hand, the energy spreading 
is much smaller than the momentum spreading and Eq. (72) or (73) is applied. Precision 
experiments of AE ~ 0, \Ap | = oo of narrow energy levels are studied with the expression 
(73), and the probability does not have pseudo-Doppler and finite-size effects then. 



4. Two-body decay: A — > B + C 

Unusual properties of the decay probability at a finite distance are studied in detail for 
two-body decay here. Decay rate is computed with S[T] and the finite-size correction to 
that computed by Fermi's golden rule is found. The correction depends on the boundary 
conditions of the experiments and is computed properly with S[T] that satisfies the boundary 
condition at T, instead of 5[oo]. 

Two-body decays of a particle A into B and C of masses m^, tub and mc satisfying 
tua > tub + mc and governed by a local Lagrangian 

J? = J? + JZ' int , (74) 
-2b = \ [{d¥A) 2 ~ m 2 A ip 2 A + (dip B ) 2 - m%ip% + (d<p c ) 2 - m 2 c <p c ], 

■&int = gtPA(x)(pB(x)¥c(x), 

in the wave zone are studied in the lowest order of coupling constant g. The characteristic 
features of decay amplitude in the wave packet scattering is seen from Fig. (2) that shows a 
space-time picture of the decay of A of a large size to a large B and C of a small size. Because 
the interaction occurs in the finite region where these waves overlap, the conservation laws 
of the kinetic energy and momentum are modified from those of plane waves. 



4-1. Average energy in the wave zone 

A kinetic energy of the wave function at a finite t is not a constant. A state vector evolves 
with the Scrodinger equation of a total Hamiltonian composed of the free part Hq and the 
interaction part H int derived by the Lagrangian, Eq. (74) 

ih^Mt)) = (H + HtofiMt)). (75) 
From perturbative expansion, a solution satisfying the boundary condition |^(0)) = \ipi) is, 

|tf(t)>= ai(i)hh) + hfe>, (76) 
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Fig. 2 The decay amplitude of A to B and C which are expressed with wave packets of 
different sizes are represented. They interact in a region where they overlap. If the size of B 
is smaller than others, the region is mainly determined by B. 



where \tpi) is one-particle state composed of A of a momentum pa and a kinetic energy Ea 
and 1^2) is a two-particle state composed of B and C 



\Mt)) = e~£ t \A,p A ), 

\Mt)) - 



* at'Eio^l^t')) = f dt '\ B ,c)(B,c\^p-\Mt'))- 







ih 







ih 



In the lowest order in g, 



ai{t) = 1. 



An energy expectation value is 

_ (g(t)jgig(fl> 

* ; " <tt(t)l*(t)> ' 

(tf (t)|£T|tf (t)) = |a a (t)| 2 ^<Vi|^i> + 2Re [ai(i)(^ 2 |tfmt|V>i>] + (£ B + £ C )(^ 2 |Y>2>, 
(^(t)|^(t)) = |a 1 (t)| 2 ^i|Vi) + (V'2|V'2). 
At the infinite t, 



(77) 
(78) 

(79) 

(80) 



{B,C 



H int (0) 

ih 



a 1 (t)(i; 2 \H in t\^i} = -(l-e^ t / h ) 
to 

uj = Ea — Eb — Ec, 



\A) 



(B,C\ 



ih 



\A) 



O(l) 



(81) 
(82) 



hence the expectation value of the interaction Hi n t(0), Eq.(82), is negligible compared to 
Eq.(81) and the kinetic energy as well as the total energy is Ea- Now for the finite t, the 
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expectation value of Hi n t(0) is not negligible. An average over a finite-time interval gives 

2 



2Aver(Re[a 1 (t)(i/j2\H int \ip 1 )]) 



{B,C\^1\A) 2 



{B ,C\^4^\A) 
in 



h 2 

Aver{(ip 2 \H Q \i> 2 )) = 2(E B + £ c ) — 



(B,C\^^-\A) 
in 



(83) 
(84) 
(85) 



The expectation value of the total energy becomes 

Aver(H) = E A - (86) 
Thus the average energy coincides with the initial energy, but the average kinetic energy is 

Ave ' iim irw^Wm — ' m 

and is different from the initial kinetic energy. Hi n t causes a transition of A to B and C, and 
is non-diagonal in the base of eigenvectors defined by Hq. Thus H% n t does not contribute to 
the total energy in the infinite t, but at a finite t, the state |\l/(t)) is a superposition of \A) 
and \B, C) and Hi n t has a finite expectation value. The total energy is always the same but 
the expectation value of Hi nt is finite in a finite t. Hence the state becomes a superposition 
of different kinetic energies and the kinetic energy is not a good quantum number in this 
region. 

E A is real in the lowest order of g and has an imaginary part in the second order, which 
represents a life time of A, ta- In t <C ta, the imaginary part of Ea is negligible. For a 
self-consistent treatment of the decay process, we start from Ea of an imaginary part and 
compute the decay amplitude and probability. The decay probability is proportional to T in 
T <C ta and becomes unity at T 3> ta- 

4-2. Transition amplitude and decay probability 

Next we study the transition probability at a finite distance. A decay of a particle A at a 
space-time position (Xa,Ta) into particles B at (Xb,Tb) and C at (Xc,Tc) in the most 
general case of the symmetric wave packets 



ga, cr B , cr c , 

of the four-dimensional momenta and masses 

(Ea,Pa] m A), (E B ,PB;rn B ), (EcPc^c), (89) 

is studied here. A life time of A expressed with an imaginary part of Ea is assumed negli- 
gible in majority places of the present paper. From the interaction Lagrangian Eq. (74), the 
transition amplitude is expressed with an integral over (t, x) 

M(A -^B + C)=g J dt J dxe-^- So)2 -^ {t - tor e R+ ^ 

= g(27ra s ) 3 H27ra t )h R+i ^9(X i ,T i ), (90) 

for finite values of a s and at, where #(Xj,Tj) denotes the conditions that to is the inside of 
the time region defined from the boundary conditions and we omit to write it hereafter. a s 
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1 




1 






+ 




+ 




a s 


OA 




OB 




o-c 


1 


V 2 


+ 


V B 
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V 2 

V C 


<*i 


OA 




ob 




OC 



and at are given in the expressions 

(91) 

4- + - + ^V- (^) 

\oa o B ocj 

The center position xo(t) is 

x (t)=x o (0) + v o (t-t o ), (93) 

of an average velocity vq, 

(V A V B V C \ . , 

vo = os [ 1 1 . (94) 

\OA OB Oc J 

R and (j) in the exponent are obtained from Eqs. (50), (51), and (53) and are given as, 

R — -^trajectory ^momentum! (95) 
^trajectory = - £ g + ^ + 2„ (v , (96 , 

^momentum = ~^{oE - VQ-5pf - y (Sp) 2 , (97) 

SE = E a {pa) ~ E b {pb) - E c {pc), $P = PA - VB - VC, 

and is a function of the momenta pj and positions Xj. 

Since trajectory is a function of the momenta and coordinates, we write it as 
^trajectory {Xa , Ta; Xi,Ti), where I stands for B or C. This is invariant under the translation, 
Eq. (52) 

-Rtrajectory(X4 + d,T A + 5] X[ + d,T [ + 5) = -^trajectory POl, ^A', Xl, T|), (98) 

and is a function of the variables Xi = Xi — v{Ti. Hence the left-hand side of the above 
equation is written as, 

-Rtrajectory(X4 + d - V A 5, T A ] X + d - V[5, Tj). (99) 

Choosing d = va$, we have the identity: 

^trajectory (Xa, ^A', Xi + Va6 - V[5, Tj) = -^trajectory (Xa, ^A] Xi,Ti), (100) 

and 



d_ 

85 

The probability is the integral 



^trajectory (X A , T A ] X l + V A S - V t 6, T t ) = 0. (101) 



\A4\ 2 does not depend on 5 from Eq. (101), and the phase space is reduced to that in 5 
component and the orthogonal components, Xt, 

J dSUdX^-^lMl 2 . (103) 
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1000 



Fig. 3 The probabilities of decays at rest Eq. (104) for the wave packet in two-body decays, 
J/V -)• M(2981) + 7 (red solid), J/^ -)• r/(1409) + 7 (green dot), vr -)• fi + v (bule dot), and 
/U — ?• e + 7 (magenta dot) . Wave packet of another daughter is 00 and that of the parent is 
cr p arent m n = 10000. The horizontal line shows the size of wave packets of the light particle in 
units of (rmj and the vertical line shows the deviations of the rates for wave packets over the 



rates for the plane waves, 1 — 
due to numerical uncertainty. 



-. Errors for 7r and \x are slightly larger than others 



The parameter 5 is not measured in the ordinary experiment and is integrated. From the 
integration over 5, we have 

Thus the probability in the system of finite a s and at is proportional to time interval, T. Its 
magnitude is independent of the parameters of wave packet from the completeness Eq. (40) 
and agrees with the value obtained with S[oo] defined with plane waves combined with ie 
prescription. 

In Fig. (3), the rates computed with the wave packets of various sizes are compared with 
those of the plane waves in various decays, J/^> — > M(2981) + 7, J — > 77(1409) + 7, 7r — )■ 
fj, + v and /i — > e + 7 which will be discussed later. Wave packet of another daughter is 00 
and that of the parent is o" paren tm^ = 10000. The value is the same for all processes. Within 
small errors, they agree. 

4-3. Various cases of wave packets 

We study the amplitude and probability of the systems (1) a s = finite, at = finite, (2) a s = 
finite, at = 00, (3) a s = 00, at = 00 in the following. 
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4-3.1. Finite a s and finite at ■ When a a, ob, and oq are finite, a s and at are also finite 
and the integrand in Eq. (90) decreases fast at t — > oo and \x\ — > oo and integrals over t and 
x converge fast, and the results of Eqs. (50), (51) and (53) are applied. The total probability 
is obtained by integrating the momentum and position of Eq. (102), and does not have a 
finite-size correction at a macroscopic T. 
When a a and as are finite and ac = oo, at and a s are finite generally. We have 

(105) 
(106) 
(107) 

The integrand of Eq. (90) decreases fast at \x — x$\ — > oo and the integral over x converges 
fast, at is finite when va ^ vb and the integrand decreases fast at \t — to I - > oo and the 
integral over t converges fast. We have 

{(Xb-Xa)^ 

Rtrajectory = —, j r , (108) 

l[pA + &B) 

(Xb - Xa)t = {Xb - X-a) - 0^ • (Xb - Xa)j=t 



1 


1 1 


1 




<?A OB 




1 


(va - vb) 


2 


Of 


oa + o B 


3 




( . 
= a s \ h 

\OA 


V B 
OB 



\V B -VA\ \VB-VA\ 

Thus the probability depends on the transversal components of coordinate Xb — Xa but 
not on the longitudinal component. The coordinate of B is integrated over the transversal 
and longitudinal components 

J dX B e m ™^y = J d (x B _ X A ) T d{X B - X A )Le 2R ^'^, (109) 

where the former variables are integrated in the form: 

J d(X B - X A ) T e 2R ™^y = it {a a + <r B ), (HO) 

and the latter variable is integrated using 0(Xj,Tj) in Eq. (90), as 

J d{X B - X a )l = \v B - v A \ j d{T B - T A ) = \vb - va\T- (Ill) 

Thus the probability is proportional to T, and does not have the finite-size correction. 
Rmomentum is expressed with Eq. (97) or with the energies of the momenta 

ob 

PA=PA : (PA-PB -PC), (H2) 

OA + Ob 

PB=PB^ — — {PA-PB ~Pc), (H3) 

OA + &B 

Pc=Pc- (114) 

Other cases of two wave packets and one plane wave are equivalent to the previous case. 
In case of vb = va] 

In the limit v*b —> va, at diverges and may cause a large diffraction effect. Nucleus trapped 
in matter have the momenta pa = Pb = and Mossbaur effect is such phenomenon that 
occurs through absorption of a gamma ray by a nucleus. 
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PC PB/ B (detected) 



F({p A -pef) 



A (parent) 



Fig. 4 Form factor in (A\J B (0)\C). 



4-3.2. Finite a s and infinite at- In finite a s and infinite at, the wave functions of initial 
and final states overlap in a long strip region, accordingly the probability shows unusual 
finite-size corrections. 
A: Small mass 

We study next the situation where the particle A and C are described by plane waves 

a a = oo, ac = oo, (115) 

of the momenta pa and pc and B is described by the wave packet of the size as and 
momentum ps- B is assumed to have a small mass mg, A is prepared at Ta and B is 
detected at a space time position (Xb,Tb). Obviously the parameters of Eq. (46) become 

a s = a Bl (116) 

Vq = V B , (117) 

1 = ^-^ = 0. (118) 

Since at = oo, the integrand in probability does not decrease with t and may receive a 
finite-size correction. 
The transition amplitude is expressed in the form, 

M = j d A xN lWB (j>B,Z)e- ipA - x+ipc - x F{{pA-pc) 2 ), (H9) 

where N\ = ig (mcms/ EqEb) 1 ^ 2 and the coefficient Nb in wb(pb,x) is Nb = (<tb/7t) 4//3 , 
F((pa — Pc) 2 ) is the form factor, and the time t is integrated over the region Ta <t< Tb- 
as is estimated using the size of a constituent object in a target which B interacts with. 
The coordinate x is integrated next and the amplitude becomes finally 

M = N 1 N B (2n/aB)H^B)^e-^ EBTB -P B ^^e-^ A -P B -P c ^ 

x [ T dte-^ EA - Ec - EB ~^ A ~ r!B ^ c ^ t F{{p A - pc) 2 ) 
Jo 

= NiNb (2n/a B ) t (47ra B ) 3 e ~ l{EBTB " pV * B > e~ ^ ®* -*0 a 

X %,-p C f)^, (120) 
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where u is 

u = E A - E c - E B - (pa ~ PB ~ Pc) ■ v B , (121) 

Because the magnitude is inversely proportional to u, M gets contributions from small and 
large u regions. The amplitude gets a large contribution at a large T from the region, 

u w 0. (122) 

A normal root of satisfying 

E A -E C -E B ^ 0, pa — Pb ~ Pc ~ 0, (123) 
and a new root of satisfying 

E A -Ec-E B ^ 0, pa-Pb-Pc^ 0, (124) 

exist. Because the kinetic energy and momentum are different from those of the initial state, 
the secondary root gives the finite-size correction due to the diffraction. Dependence of the 

did 

dp E 



amplitude on p B is determined by a root of u = and its slope Jnr-. 



Assuming \p A — Pc ~ Pb\ is small, we have 

E B (p B ) + (pa - Pb -Pc)-v B = E B (p A -pc), (125) 

and 

u = E A (p A ) - E c (pc) ~ E B (p A - pel (126) 
The probability integrated over pc becomes 

NlNl(2,/a B f(^B) J ^e-^-^r (^^T) * F{{pA _ pc?) 



sin(wT/2) 



2 



= NfN B (2Tr/a B y(4ira B ) J du { \ ) ( 12T ) 
where the spectrum density p(uj) is 

pf U ) = [ ^fce-vaWA-pB-po)* ( 128 ) 
J (27T) d 

x 5{u - E A (p A ) + E c (p c ) + E B (p A - Pc))F((p A - Pc?)- 

Because p(0) is finite and p(ui) decreases rapidly at the large u region, as is given in Appendix 
A, the following integral converges at a finite T, 

sin(o;T/2)\ 2 , x „ f„ ,„ N 1 



<M u ) PM = T|2^(0) + -C|, (129) 

where £ is equal to C(T) in Appendix A. Thus the finite-size correction becomes finite. 

The finite-size correction to the total probability integrated over whole momentum region 
of pc is obtained easily with the correlation function [21], 



J (27rW' ' = ~E~b J A A ,c(5x)e 



i<j>(Sx) 



(130) 



where N 2 = g 2 (4ir/a B f /2 V~\V is a normalization volume for the initial state A, x? = 
Xb + v B (ti — T B ), 5x = x\ — X2, 4>(5x) = p B -5x and 

A^(fe) = (^3 / ^f-^-^niPA - Pc) 2 )- (131) 

In the right-hand side of Eq. (131), the integration region of the momentum pb is that of 
the complete set and is reduced to the smaller one if the integrand |A4| 2 vanishes in some 
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kinematical region. This happens for the amplitude of plane waves at the asymptotic region 
T = oo, which includes the delta function, 5^(5p), from the integration over x reflecting 
the conservation law of kinetic energy and momentum. The phase space of the final state 
becomes proportional to the initial energy, then. Now in the right-hand side of Eq. (131), the 
coordinates are fixed and are not integrated. Thus the correlation function Aa,c($ x ) does 
not include 5^(5p), and p c is integrated over the whole region.. 

Because the probability is finite, integration variables can be interchanged. For fh? = 
m 2 A — mfj > and a real pa, [18, 19, 21], and from Appendix C, 



Aa,c(&c) = 2i 



F(-m 2 )D A ( -i f^J(A) + }, h „ t ] 



d8x J \ Air 
d \ ^ 1 / . d \ d 



(132) 



d5x I ' 11 V V ddx J dm 2 

• ~ 2 ■ ~ 9 

f snort = —ei-Xiimd) - fc(«)Ji(e)} - £^° ix)Kl{0 > 

where e(6t) is equal to = +1 or —1 for positive or negative St, respectively, A = (5x) 2 , 
£ = m^f\, and Ni, J\, and K\ are Bessel functions, j short na s a singularity of the form 1/A 
around A = and decrease as e _m V / l^l or oscillates as e' m V^I at large |A|. The condition for 
the convergence of the series will be studied later. The formula for A with a finite life-time 
is obtained later. The last term is 

1 f nv^ _ \2w_ o n \ d PC 



F((pA-pc) 2 )9(PA-Pc)^Te-^-^> 5x . (133) 



(2vr)3 7 — ^ ^^(pb) 



For m 2 = m 2 ^ — m 2 -, < 0, 



Aa,c(&c) = 0. (134) 



Thus Aa : c($ x ) is composed of the light-cone singularity 5(\)e(5t) [18, 32], the regular 
terms written by Bessel functions, and li- The former two terms come from the integration 
from Ea < Eq, and are finite in a finite T. So using this expression the finite T correction 
which is unable to obtain with standard calculations of plane waves can be found. Because 
the integration region for this is outside of the kinematical region of conserving the energy 
and momentum, this integral vanishes at T = oo in fact. I2, on the other hand, comes from 
the region Eq < Ea, which is the kinematical region satisfying the energy and momentum 
conservation, determines the quantities at T = 00. This expression of writing the probability 
with the light-cone singularity converges and is valid in the kinematical region 2pa-Pb < ™c-> 
where vr? c = m 2 A — vr? c . 

Substituting the expression of Aa : c(Sx) into Eq. (130) and integrating over x\ and X2, we 
have 
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for the leading singular part and 

J 1/A = J d£ 1 d£*W^e-^(* 1 -* B -** tl - 1! ' ) y 



4vr 2 A 



a B ( \ 2 _- ffB |fi,|» 1 J$ c (St) 



for the next term of the form 1/A. 

Finally we integrate t\ and t 2 over the finite region T = Tb — T^, and we have the slowly- 
decreasing term g(co B T) 

i [ T dt 1( it 2 e -^e^ 5t = T(g(cj B T)-7v), (137) 
Jo 1**1 

2 

L0 B = E B - \PB\ - 



2E B ' 

and the normal term Go- <?(cj_bT) is generated from the light-cone singularity and related 
term and satisfies 

5(0) = 7T, (138) 
s(w B T) -> — T^oo, (139) 

and vanishes at T = oo. Go is from the rest 

Go = 2^ | ^|t*(£a - £b - J Bc(pc))e- CTfl( ^- fe -^ )2 0(^ - (140) 
and conserves the kinetic energy and momentum approximately 

PA- PC = PB, (141) 

and gives the asymptotic value. Due to the rapid oscillation in St, Go gets contribution from 
the microscopic \St\ region and is constant in T. Integration of this term does not depend on 
a B and agrees with the normal probability obtained with the standard method of using plane 
waves. In the region 2pA-p B > fhc, ^-A,c{dx) does not have a light-cone singularity and the 
diffraction term exists only in the kinematical region 2pA'PB < ™c- This region depends 
upon the mass of G, hence the diffraction terms of all three mass eigenstates converge in the 
union of the kinematical regions of the three masses, 2pA-pB < ™c- ^ ne diffraction terms 
exists in this region. 
We have 

/ ^\M\ 2 = §^{F(-rh 2 )g(u> B T)+F(m 2 B )G }, (142) 

where N 2 = g 2 (4ir/a B ) 3/2 V^ 1 , V is a normalization volume for the initial state A. Form 
factor gives different corrections to the diffraction and normal terms. They are evaluated 
later. 

B: massless particle m B = 

For a massless B, the leading singularity 5(X)e(5t) cancels on integrating over the times 
and the next term proportional to 1/A gives a dominant contribution. The integral of this 
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term 



leads 



J 1A = / d ft( JS 2e i *( fa ) e -^a( £ '--*-- t; -<''- T ->) 3 i -L-, (143) 



2 V^bIP-bI 2 / 1^1 

This term also has the universal dependence on \St\ and its integration over the times 
becomes, 

f dhdt 2 J 1/x = (a B 7r)^ ( 1—) 2 e — l^l 2 f dtl dt 2 -^—. (145) 

J 2 \7ra B \p B \ 2 J J \h-t 2 \ 

The integration over the times in a finite region from e to T is 

J dhdh * = T (^2 log ^ - 1^ , (146) 

and 

= <«*> ! T (dfeF) 5 e ~"'"'° T ( 21og 7 " ' (147) 

This term gives the probability 

P dif f = N 3 Jdp B e-°^ 2 , (148) 

where 

A 3 = 8T 5 2 f^) f 2 log--lV (149) 



4 7 V e 



Large time: T > ta 

If T is larger than the life time of A, ta, Eqs. (138) and (146) are replaced with 

f' T e(5t) p ti+«a 

i/ dhdh^-fe^e -a =g{u B ,T;T A )-g(uJ B ,T;TA), (150) 



g T (u) B , oo] t a ) = 0, u B = E B - \p B \ 



m 



2E B ' 
and 



1 t 1 +t 2 

dt 1 dt 2 -, re t a . (151) 

\h — 1 2 \ 



N3 becomes approximately 



iV 3 = 8r Aff 2 ^(21og^-l). (152) 

Thus, the system of at = 00 has the finite-size correction of the form Eq. (142) for T <C ta, 
and g(Tuj B ) in Eq. (142) is replaced with g(T,w B : ta) in T w r^. The correction depends 
on T in the universal manner and on the size of wave packet in magnitude. At a B = 00, the 
correction becomes infinite. 
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We compute the total probability next. From the integration over Xg, the total volume V 
is obtained and cancelled with the normalization of A. The total probability, then, becomes 
the integral of a sum of Go and 3(wbT), 

^3 / [F(-rh 2 )~g(co B T) + F(m 2 B )G ] , for T « r A , 

J ^ (153) 
N 3 J JpJ^B [F{-rh 2 )~g{T^ B -T A ) + F(m%)G ] , for r A < T, 

where N3 = 8Tg 2 o~B- The second terms, Pnormah i n the right-hand sides of the Eq. (153) are 
independent of T and o~b and agree with the standard value computed with the plane waves. 
<7(wbT) and <?(T,cjb; t a ) in the first terms depend on ujb and T, and are corrections due to 
the finite distance between the initial and final states. Their magnitudes of the first terms, 
Pdif fraction, at T -> oo, are proportional to 

Pdif fraction = NF(-m 2 )-^ = NF(-rh 2 )^^, (154) 

where N is constant. Pdif fraction becomes significant for large {o~bEb) /m 2 B , i.e., small mass 
or large wave packet. 

4-3.3. Infinite o~ s and infinite at- When three particles are plane waves, a s = at = 00, the 
scattering amplitude and cross section are the standard one if i/; nt (£)e- e l*l is used. The space- 
time coordinates (t, x) are integrated over the whole region and the energy and momentum 
are strictly conserved. The asymptotic values thus obtained with 5 [00] 

M = (2vr) 4 9 5^ {p A -PB- p c )f, (155) 
P = 9 2 \f\ 2 x {phase space), (156) 

agree with the asymptotic values obtained with S[T]. If the convergence factor e _e '*' is 
absent, the limit T — > 00 is not unique and is consistent with the diverging correction in 
ob — > 00 of the previous case. 

4.3.4. Coherence length. The coherence length found from the amplitude of the initial and 
final states expressed with wave packets is finite. From Eq. (91), the integral in x converges 
for a finite a s and that in t converges for a finite at- at becomes infinite with v A = vb and 
ac = 00, or o A = ob = 00. In the latter case, the coherence length is KEc /(m^c 3 ). 

4-3.5. Asymmetric wave packets. For the asymmetric wave packets, the integral over (t, x) 
is expressed by 

where the sizes of the Gaussian exponents and other parameters are given by complicated 
expressions. Experiments of 5E <C \5p\ are studied with the asymmetric wave packets. 
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5. Emission and absorption of light 

Radiative transitions of particles 

A + ^^C, (158) 

expressed with wave packets are studied in various parameter regions. Electromagnetic 
interaction is expressed with 

H^t = e J dxJ n (x)A^(x), (159) 

where A n (x) is a photon field and Ju(x) is a electromagnetic current. A matrix element of 
the current between eigen states of energy and momentum are written as 

(C;p c \Ju(x)\A;p A ) = e^-P^ {C; Pc \Ju(0)\A;pa) , (160) 

where 

(C; Pc \M0)\A;pa) = r u F((p A -p c ) 2 ), (161) 

with the form factor F((pa — Pc) 2 ) and the spin dependent factor, T n . We assume one form 
factor for simplicity, but it is straightforward to extend to a case of many form factors. In 
the normal term of the radiative transition, the energy-momentum is conserved, and 

F((p A -p c ) 2 ) = F(k 2 1 ) = F(0), (162) 

hence the coupling strength is determined by ^(0). 

In detectors, fundamental processes of a photon are, photo-electric effect, Compton effect, 
or e + e~ pair production. Wave packet sizes of the photon, <r 7 , are nuclear sizes for the pair 
production due to nuclear electric field or atomic sizes or larger for the photo-electric effect 
and Compton effect, depending on the energy. 

5.1. Universal background 

The transition probabilities of radiative processes receive the finite-size corrections under 
certain situations and their energy spectra are modified by pseudo-Doppler effects. Since 
the finite-size correction is caused by the states that violate the conservation law of kinetic 
energy and momentum, the corresponding events look like backgrounds even though they 
are produced dynamically. They have universal properties and magnitudes that depend on 
the experimental apparatus. 

5.1.1. Universal background. The universal background derived from the finite-size 
correction resulted from 

e^ T - 1 



± 0, (163) 



UJ 

is inevitable consequence of Schodinger equation. Since they are generated by the states 
of the kinetic energy different from that of the initial state, that is positive semi-definite 
from Eq. (35) and is added to the normal component in the wave zone. Its magnitude is 
computed rigorously in relativistic systems Eq. (153). The correction vanishes in the particle 
zone. Energy spectrum for wave packets is distorted in both particle and wave zones due to 
the pseudo-Doppler effects even though the total probability agrees with the normal value. 
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5.1.2. Form factor . Nucleus, atom, and molecule are composite states and have internal 
structures. So they have finite extensions and interact with photon or neutrino non-locally. 
This non- locality is negligible if the size R and the photon momentum fc 7 satisfies k^R <C 1, 
where the multi-pole expansions are applicable. 

For X-rays from atoms, they are about 

k y R = 10~ 3 ; A; 7 ~ keV, R = KT 11 [m], (164) 
and for transitions of nucleus 

k^R = KT 1 ; A; 7 ~ MeV, R = 10~ 15 [m]. (165) 

Since k-yR is small, 

F(( PA -pc) 2 ) = F(0). (166) 

5.1.3. Life time effect. If the parent A has a finite life time, ta, that modifies the results. 
In a region 

CT A < y/o~A, (167) 

the integral over the times in the transition probability gets a dominant contribution from 
the region 

t < ta- (168) 

Then the effect of the wave packet is diminished and the pseudo-Doppler effect becomes 
negligible. If the life time satisfies 

cta > V(ta, (169) 

the integral over the times in the transition probability gets a dominant contribution from 
the region 



t < (170) 

c 

and the pseudo-Doppler effect is prominent. 

5.1. 4- Photon effective mass. Photon is massless in vacuum but its property is modified 
in matter due to dielectric constant. In high-energy regions, the refraction constant behaves 
with the frequency as 

ra = l--|, ( 171 ) 



where u p is the plasma frequency and is given as 

o NZe 2 , . 

<4 = — > 172 

y e ?n e 

depends on material, density, and other parameters. The wave vector satisfies 

{ckf = u? - W p, (173) 

and the energy dispersion becomes 

p = E{p) 2 - (hio p ) 2 . (174) 



30/57 



Thus the photon has an effective mass 



NZe 2 , . 

m eff = fu — — , (175) 

V e m e 

where N and Z are a number density and atomic number of gas, and m e is the electron's 
mass. The m e ff depends upon the density of matter and varies. The high-energy photon 
behaves like a massive particle. 

5.1.5. Light-cone singularity for general systems. For particles A and C of internal struc- 
tures, Eq. (161) is substituted to Eq. (131). As is shown in Appendix C, the singular part of 
correlation function is written in the form, 

A fx \liqht-cone t?( _2 , 2 \ a (0), light- cone , x 

Aa,c{ox) y = F(—m A + m c )A AC (dx), (176) 

where A^ c (5x) is that of the point particle. Thus the form factor 

F(m 2 c -m A ), (177) 
determines the strength of the singularity and are given in Appendix C, as 

0(1); hadron, positoronium, light nucleus, 
O(10 _1 ); [iN atom, heavy nucleus, 



F(m" c - m A )/F(0) 



(178) 

0(10 5 ); /ie, K — electron, 
O(10~ 10 ); atom, 

Thus the form factors do not modify the magnitude of light-cone singularity for hadrons, 
light nucleus, positoronium, and reduce to 1/10 for //-N atom and heavy molecules. For 
/i-e, K-electron, and atoms, the magnitudes become extremely small. Eq. (177) is almost the 
same as the on-shell coupling, Eq. (162) in the former but much smaller in the latter. Because 
the singularity is caused by the waves of translational motion which retain the relativistic 
invariance even for the particles with internal structure, but the magnitude depends on their 
sizes. 

5.2. Emission of lights 

5.2.1. Decay inflight in vacuum. 1. Finite o~ A and cr 7 : pseudo-Doppler effect 

The amplitude of a radiative decay of A to C and a photon 7 of momenta, positions and 
wave packet sizes 

A : (X A ,E A ,p A ,a A ), 

7 : (Xy, Ey,p y , ct 7 ), E 2 -p^ = 0, 

C:(pc,E c ,o-c = 00), (179) 
is expressed with the matrix element of the current operator and the photon field 

M = j d 4 x(O|J At (x)|A)( 7 |^(x)|0) 

^ Xe i(PA-PC~P-,)-X p^^Py-Xj-^-iX-Xy-Vyit-Ty)) 2 

X e -iPA-X A -^(x-X A -v A {t-T A )f _ e R+i<t>^ / 18Q x 



Fab = (C\J^(0)\A)e»(p 



31/57 



where e^(A; 7 ) is a polarization vector of the photon. We have \Ai\ 2 in the form 
\M\ 2 =N 2 f ^^^e 4 ^-^-^)-^-^)-^^^-^-^^-^)) 2 



4 „J 



x e -^rS(*-^(^Tx))v y(PiliJJo)( ^ _ Ml), (i8i) 

p 7 

in Coulomb gauge, 

A(x) = 0, V -A{x) = 0, (182) 
where N is the normalization factor, and 

W itj (p A ,pc) = {ClJiMlAWClJjWlA))', (183) 

i j 

{S iJ _P^L ) = J2e%)(e^k J )T. (184) 

P-y 

R in Eq. (180) is composed of the momentum-dependent part Rmomentum and the 
coordinate-dependent part Rtrajectory The former is 

Rmomentum = -^(Ea(pa) ~ E C {pc) ~ E^^)) 2 - ^{PA ~ PC ~ Pjf , 
Pj=p 7 +—(PA-Pc-Py), (185) 



where a s , at and vq are 



*. = (186) 



a A + cr. 



7 



wo = — ^7 — w 7 , (187) 



a a + cr, 



7 

I = 5_3?=_£_. (l 88 ) 
cr f cr 7 cr s cr^ + cr^ 

Thus, the energy-momentum satisfies the modified conservation law. The momentum is 
conserved approximately around the center 5p = 0, whereas the photon's energy at the 
momentum p 7 fulfills the approximate conservation law. Implications will be studied in 
detail shortly. 

The position-dependent exponent is written in the form: 

X^ I XA Xry \ 

Rtrajectory = ~ ^ " ^ + 2<7 S ^— + _J 
/ =* \ 2 



2(7,4 2(7. 



1 



1 



(X A - X 7 f - ^(i/ 7 • (Xa - X^f 



2 _^„-? ^ v . v ^^2 

"7 



2(17,4 + OVy) 



(189) 
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The probability is expressed as 



P — N exp {^Rmomentum ~\~ ^Rtraj ectory) I Ea,B \ 



N 2 \F A , B \ 2 exp 



-a t (E A - E 7 (p 7 ) - E c ) 2 - cr s (p A -p c - p 7 f 



x exp 



(J A + 



(190) 



and has no finite-size correction. Thus the total probability agrees to that of plane waves. 
Nevertheless, the energy spectrum of Eq. (190) is distorted due to pseudo-Doppler effect. The 
photon's momentum is distributed around a center pa — Pc an d the photon's energy at a 
momentum p 7 is distributed around Ea — Ec- If <J S is small and at is large, the momentum 
distribution is wide but the energy E(p v ) coincides almost with Ea — Ec- The observed 
photon's energy is E^(p^) and is given from Eq. (185) 



E i(Pi) = E i(P 1 



(j - 



(PA-PC-P-I)) 



a. 



Ea — Ec # 7 • (pa - Pc - Pj)- 



a 



(191) 



Thus _E 7 (p 7 ) is very different from Ea — Ec- 
Photon is on mass shell and satisfies 



E&f-p^O. 



(192) 



In an event where the energy-momenta (Ea,Pa), (Ec-,Pc)i an d (-E 7 ,p 7 ) are measured, and 
the momenta satisfy 



P-y t^PA-PC, 



(193) 



the photon's energy at the momentum p 7 satisfies 



E A - E c = £ 7 (p 7 ). 



(194) 



Consequently the mass shell condition at p 7 



E 7 (^) 2 - ~p 2 = (E A - E c ) 2 - ~p 2 = 0, 



(195) 



is satisfied. Substituting jj 7 , we have 



(E A - E c ) 2 - (p 7 + —(pa-Pc- P 7 )) 2 = 0, 



(196) 



which gives a relation of the energies and momenta with the ratio a s /a^. Measuring the 
energies and momenta, the ratio a s /a^ will be determined. 
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In a situation of 

cr 7 < cr A , (197) 

we have 

<7 S = cr 7 , crj = <7a, cr s <C <7t, (198) 
v = — — p 7 = — ^ — p-y + (pk - pb). (199) 

The central values of energies and momenta satisfy 

(Eyfa)) = (E A -Ec), (200) 
(Py) = iPA-pc), (201) 

with variations 

8E = (202) 



\5p\ = -L, (203) 



The energy spreading is narrower than the momentum spreading, 

SE < |<5p|, (204) 

hence the constraint to the energy is more stringent than that of the momentum. 
Heavy A and C (Pseudo-Doppler effect combined with Mdssbaur effect) 

If C and A are a ground state and an excited state of a heavy atom which bound together 
to become massive objects, the correlation function of Eq. (183) does not vanish only at the 
same momenta, 

PA=PC, (205) 

like those of Mossbaur effect. We study the photon's energy spectrum when this condition 
is satisfied in a large wave packet a a = gc- The reduced momentum becomes 

p\ = ^ P 7 > ( 206 ) 



from Eq. (198). The energy of the massless particle is proportional to the momentum and 

(E 7 (p\)) = ^2_<£ 7 (p 7 )), (207) 

OA ~r o 7 

Substituting Eq. (200), we have the expectation value of E 1 

{ E i{Pi)) = KAE dectron , (208) 



0~A + C7 



AE clectr on = E A -E C , (209) 



which is much larger than the energy difference E A — Eq- Thus the product of average 
energy with the time interval for the light is equal to that for the atom 

a 7 (£ 7 (p 7 )) = {a A + <x 7 )(£a - E c ). (210) 

Now, £J 7 is the size of particle which the photon interacts with and a a is that of the atom, 
they are proportional to the average-time intervals of their reactions. Thus the conservation 
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law Eq. (62) for the energy is satisfied for average value. This unusual phenomenon occurs 
because the electromagnetic interaction takes place in the narrow space-time region where 
the wave functions of the A, C and photon overlap. When <t 7 is much smaller than cta, the 
region has the £1X621 (J <y and moves with the velocity u 7 furthermore. Hence the energy is 
conserved in this moving frame where the photon has the effective energy E^(p^), which is 
much smaller than E^{p^). Hence the average energy of 7 becomes much larger than the 
energy difference between A and C. This is the pseudo-Doppler effect caused by the wave 
packets. 

The condition, Eq. (197), is fulfilled in various situations. A molecule in gas propagates 
almost freely and an atom is bound in solid. The wave packet size of a molecule in gas is 
given by the square of mean free path and is of order of 10 -14 [m 2 ], whereas that is the 
atomic distance in a solid of the order of 10 -20 [m 2 ]. Hence we have 

k = — = 10 6 . (211) 

Consequently the photon in this situation interacts with the atom in solid with the energy 
KAE e i ectrori . For AE c i ectron = 0.1 [eV], Ey can be as large as 100 [keV]. Various anomalous 
luminescence of X-rays or 7-rays [33-35] may be connected with this energy enhancement. 

For a photon produced from excited atom in solid and interacts with a nucleus in a solid, 
we have 10 -20 [m 2 ] for the former size and 10 -28 [m 2 ] for the latter size, and 

k= — «O(10 8 ). (212) 

Consequently the photon produced from excited atoms interacts with a nucleus with much 
larger energy than the energy difference AE e i ectron = Ea — Eq- Because the photon-Nucleus 
cross-section is much smaller than that of the photon-atom scattering, the probability of this 
event is extremely small. 

A similar phenomenon is expected when charged particles propagate in magnetic field. A 
plane wave of the charge q and mass M in the magnetic field B 

e -i(E(p a ) + ^-p )(t-T )+ipa-x ^ ( 2 13) 

has the phase proportional to Cyclotron frequency 
ihave like plane waves in a ti 

and proton are 



These waves behave like plane waves in a time region less than Tj = -jf. Tj for the electron 



2ir Mi , . 

T l = — = -^i = e,p. (215) 
q\B\ 



Thus the waves have the different sizes and its ratio is 

T R m, 1 



T p m p 2000 



(216) 



Thus, the photon emitted from the atom interacts with the electron in a magnetic field can 
reveal the same energy enhancement. 

The anomalous enhancement of the photon's energy results from the overlap of wave 
functions of different sizes. They occur when the photon's wave packets, which are the sizes 
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of the wave functions with which the photons interact, are much smaller than the parent's 
wave functions. Hence, the rate of the events may be quite low. 
2. Infinite a a and finite cr 7 :finite-size correction 

The amplitude of a radiative decay of A to C of plane waves and a 7 of momenta, positions 
and wave packet sizes 

A : (X a ,E a ,pa,cta = 00), 

7 : (X y ,E y ,p^,a y ), E 2 -p 2 = 

C:(pc,E c ,a c = 00), (217) 
is expressed with the matrix element of the current operator and the photon field 

M=eJ d A x{C\ J M (x)|A)( 7 |^(x)|0) 

=e [ d 4 xe i ( p - 4 - pc -^> x (C|J M (0)|A)e^(fe 7 )e^'^"^ (a? " XT "^ { *" T "' ))2 



(218) 



We have \M\ 2 in the form 

\M\ 2 = N 2 j d 4 x 1 d 4 x 2 e 



i 3 

xW itj (p A ,PcW' j -^). (219) 
P 2 

Integrating over pc with a variable r = pa — pc, we have 
dpc 



{2nfE c 
74, 



A PA - Pc y( Xl - X2 ) Wtj{pA ^ pc) 



= [ *ljlm J_ Wy(p Aj p A -r)e*-^-»), (220) 

J {liry \_r z — IpA • r + m A — m c — te J 

which has the light-cone singularity 

— 5(A)e(ti - t 2 )Wij(p A ,PA - Olr 2 =m=,-m| 5 ( 221 ) 

from the integration over the momentum r = (r°, r ) of the region 

r 2 - f 2 = m 2 c - m 2 A < 0, r° < 0. (222) 

It is noted that \{pa — Pc) 2 \ = l^c* ~~ m A\ ^ s sman an d \Wij(pA,PA — r)\ is almost the same 
as the on-shell matrix element of the radiative transition. Eq. (220) has also regular terms 
and one of them is generated from the above kinematical region and the other terms are from 
the region < r° < p A . The latter coincides with the normal term of the decay probability. 
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Thus we have 

d 3 p c 



/ \AA\ — Pnormal ~\~ Pdif fraction? (223) 
Ec 

Pdif fraction = Cg^LOyT), U!j = , 

where C is determined by the wave packet size. From the convergence condition in the 
expansion Eq. (220), the light-cone singularity exists in the momentum region 

2pa ■ q < m\ - m 2 c . (224) 



5.2.2. Decay at rest in solid . A decay of A in solid to C and a photon 7, which have 
momenta, positions and wave packet sizes 

A : (X A ,E A ,PA = 0,a A ), 
7 : (X 11 E 1) p 1 ,a 1 = 00), 

C :(X A ,p c = 0,E c ,a A ), (225) 
is a kinematical region of Mossbaur effect. The amplitude 

M{A^C + ~f)=g j d 4 xw A (p A = 0,x)w c (p C = 0,x)e ip ->' x , (226) 

where 



WA = N A [^y e -^-^ 2 -^, (227) 
I^Ve-^-^ 2 -^, (228) 



o A ) 

x\ = X A , (f>£ = m A (t - T A ), x% = Xa, 4>q = m c (t - T c ), 

is given as 

M (A -> C + 7) =Ne im * TA - m cTc £ dte-^ mA - mc - E -^ J d^e'^ 3 '^^ 3 

=Ne i* 2 M(™A -m c - £ 7 )T/2] ^ (22g) 

(m A - m c - £7) 

In the above equation N and <I>o are constants. The square of the modulus of T is expressed 
in the form, 

r-T 



f dtl dt 2 ^e-^ mA - mc - E ^-^e~^, (230) 
Jo E 1 



where the integral 



J_ e -K m A-mc-E~ l ){t 1 -t 2 ) e -^f-p^ (231) 



is a smooth and short-range function of t\ —ti- Hence, the total probability is proportional 
to T and has no finite-size correction. 

Particles in liquid are descried also with wave packets and the probabilities of their 
reactions are studied in the same way. 
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5.2.3. Decay in flight in dilute gas . A photon has the effective mass in X-ray of 7-ray 
region in dilute gas and the rate is modified by the large finite-size correction. A radiative 
decay of A in flight in gas to C in flight and a photon 7, which have momenta, positions 
and wave packet sizes 

A : (E A ,p A ,a A = 00), 
7 : (Z 7 ,E 7 ,p 7 ,(T 7 ), 

C: (E c ,fc,a c = 00), (232) 

is studied in a similar manner. Since a a = o"c = 00 1 the amplitude is expressed in the form of 
Eq. (218) with the effective mass of high-energy photon in X-ray or 7-ray regions, Eq. (175). 
The probability to detect this photon is given in Eq. (153) of the finite-size correction. The 
frequency that determines the finite-size correction for this photon of the energy is, 

miff 

U = Mlf ' (233) 

which gives a macroscopic distance. 



5.3. Absorption 

An absorption of 7 is studied in a similar manner as the decay process. Process of A, C, and 
7 of the parameters, 

A : (Xa,Ea,PA,cta), 
7 : (A ? 7 ,- J E 7 ,-p 7 ,c7 7 ) 

C :(pc,E c ,ac), (234) 

is described by replacing the sign of the photon's momentum in the previous amplitudes, 
Eq. (180) or Eq. (218). The distribution function deviates and central value of photon's energy 
E 1 (p 1 ) becomes different from E A — Eq by the pseudo-Doppler effect, and the probability 
receives the large finite-size corrections in certain parameter regions. 



6. Implications in particle decays 

Implications of the probabilities modified by the finite-size correction or pseudo-Doppler 
effect are studied in decay experiments. The former correction depends on the mass, energy, 
life-time, and time-interval in the universal manner and its magnitude depends on the wave 
packet sizes and the internal structures. The effect of internal wave function to the light- 
cone singularity is analyzed in Appendix C and it is shown that for hadrons, nucleus, and 
positoronium, the internal wave function does not modify the magnitude, but for atoms 
that is different. If the initial wave packets are small or T s> t, the overlap of wave functions 
becomes negligible, whereas that becomes large if the initial wave packets are large and 
T < r. Especially the probability reveals various unusual behaviors for (JinUiai S> & final- The 
finite-size effect is easily observed directly with the measurements made with the detector 
located at various L. Conversely, the latter correction becomes large for the small wave 
packets, and the energy spectrum modified due to the pseudo-Doppler effect is easily observed 
with the detector if the energy resolutions and other properties of the detector are understood 
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Fig. 5 The energy spectrum of gamma in J/^ decay at rest to M(2981) and 7 at T = 00 is 
shown. The horizontal line shows the energy of 7 in MeV for a^m^ = 14.6, 100, and 00, and 
the vertical line shows the probability. Wave packets of another daughter and parent are 00 
and dparentm^ = 10000. The probability is the same for wide region of parent's wave packet. 
The spectrum is sharp for the plane wave and broad for the wave packets. The position of 
the peak shifts for the small wave packet due to pseudo-Doppler effect. 

well. If they are unknown, the parameters of the detector are determined by comparing the 
theoretical values with the experimental data obtained from a standard sample. Calibration 
of the measuring apparatus may be used for this purpose. 

We study various decay processes and present magnitudes of the finite-size and pseudo- 
Doppler effects for the parents of plane waves and the detecting particles of wave packets. 
Life-times of the parents are included, and spin independent components are studied. 

6.1. Pseudo-Doppler effect 

The energy spectrum is modified by the pseudo-Doppler effect in wide area and the distortion 
must be known not only for a precise analysis of experimental data but also for understanding 
physical phenomena. The comparison of the rates computed for plane waves, and wave 
packets of various sizes is given for J/^S — )■ M(2981) + 7 in Fig. (5). The total rates integrated 
over the final states agree but the energy spectra differ depending upon the wave packet size. 
The distributions and the shifts become wider and larger in smaller wave packets. 

The broadening and shift of energies in other process such as J/^> — > M(2981) + 7, J/^> — > 
77(1409) + 7, 7r — > [i + and /1 — > e + 7 are compared. They are sensitive to the wave packet 
size as in Figs. (6)-(8). Fig. (6) shows the variance of the final energy of various process. They 
are almost on one line. Hence, the wave packet size can be known from the variance of the 
energy of the final state. In Fig. (7), the normalized variance of energies of final states are 
presented. Those of heavy particle are different from those of light particle. In Fig. 8, the 
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size of wavepacket 



Fig. 6 The variance of final energy, AE = y {(Eg + Eq) 2 ) — {Eb + Eq) 2 in particle 
decays at rest measured at T = oo for J/* ->■ M(2981) + 7 (red solid), J/* ->■ 77(1409) + 
7 (green dot), tt —>•// + v (blue dot), and /i — >■ e + 7 (magenta dot). The horizontal line 
shows the size of wave packets in units of am\ and the vertical line shows the variance 
, AE. Wave packets of another daughter and parent are 00 and (T paren tm 2 = 10000. Curves 
are almost on one line. 



average energies of the final states are compared with the initial energies. The deviations 
are clearly seen, and the total energies of the final states become larger than those of the 
initial states. 



6.1.1. Radiative transition of atom and positronium. An atom is a bound state of a 
nucleus and electrons and is heavy. Radiative transitions of atom, from an excited state to 
a lower energy-state, which emits a photon, are examples of the two body decays. Electrons 
bound to a nucleus have sizes of about 10 -10 [m] and energies of about 10 [eV] or less. 
The photon is detected through its interaction with matter in a detector. Among various 
reactions, photo-electronic effect is the most important reaction, where an electron is emitted 
from the interacting photon with electrons. We assume here that the electron which the 
photon interacts with is a bound electron in atom at rest. The size of its wave function is 
about 10~ 10 [m]. So cr 7 has this size. For the initial particle A, a a has either about the same 
size, (1) 10 -10 [m], for A in matter, or (2) larger than 10 -10 [m], for A in vacuum or in dilute 
gas. In exceptional situation, (3) a a is smaller than 10 -10 [m]. In experiments of 5E ~ \Sp\ 
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Fig. 7 The variance of energy of the final state over the average energy of the 
final state, AE/(Eb+Ec) in particle decays at rest measured at T = oo for J/ty — > 
M(2981) + 7 (red solid), J/V -)• 7/(1409) + 7 (green dot), it -)■ fi + v (blue dot) and /_i -)• 
e + 7 (magenta dot). The horizontal line shows the size of wave packets in units of erm 2 
and the vertical line shows the deviations of the energies. Wave packets of another daughter 
and parent are 00 and (y varent m^ = 10000. AE/(E2 + E3) is proportional to (am 2 ) -1 / 2 . 

of the following three cases of the wave packet sizes, 

(1) : a a ~ c 7 , 

(2) : a a > o- 7 , 

(3) : a a < o- 7 , 



the energy spectra are modified differently. 

Positoronium is a bound state of an electron and its anti-particle, positron. Positoronium 
of positive charge conjugation decays to two gammas and that of negative charge conjugation 
decays to three gammas. The former is a second order QED process and the latter is a third 
order QED process, and the phase space are also different. Hence their decay rates are very 
different. 

6.1.2. J/ty radiative decay. Photons produced in the decay 

J/m M + 7, (235) 

have energies in GeV region and may receive the pseudo-Doppler effect. The J/ty is produced 
in e + e~ reaction and has a size determined with beam sizes, and meson M is detected with 
its decay products, which are stable hadrons such as pion, kaon and others. These charged 
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Fig. 8 The deviation of average energy of the final state from the initial 
mass, m,A — {EB + Ec) [MeV], due to the pseudo-Doppler effect is given. There 
are finite deviations in various decays at rest measured at T = oo for J/^f — > 
M(2981) +7 (red solid), J/^f ->• r/(1409) + 7 (green dot), vr ->• \i + v (blue dot), fj, ->■ e + 
7 (magenta dot) and ip' — > M(2981) + 7 (light-blue dot). The horizontal line shows the size 
of wave packets in units of arn^ and the vertical line shows the deviations of the energies. 
Wave packets of another daughter and parent are 00 and a par ent'm'i = 10000. 

particles have semi-microscopic sizes and o~m has the same size. 07 is the order of the nuclear 
sizes. 

These processes are important for QCD dynamics for M = cc state [36] or for glueball 
M = glue ball [39] . Magnitudes of the corrections to the probabilities are not negligible as 
are presented in Figs. (6)-(8). The decay 



is almost equivalent to Eq. (235) except the phase space and has smaller pseudo-Doppler 
effect due to large 7 energy. Experiments show a difference between Eqs. (236) and (235) 



6.1.3. Two gamma decays of a heavy scalar particles. Positoronium, neutral pion, charmo- 
nium P-states, and Higgs scalar decay to two photons. They are identified with reconstructed 
photon's energies and momenta. Detection of photons are made with photo-electric or Comp- 
ton effects in low energy regions and with the e~e + pair productions in high energies. The 
bound electrons of atoms in insulator have the size of 10 -10 [m] so 07 for the former processes 
are of this size. The e~e + pair is produced by an electric field around a heavy nucleus, which 
is short-range of nuclear size. Hence the wave packet size for the latter process is approxi- 
mately the nuclear size in the high energy. Hence, the wave packets sizes of the gamma vary 



ij/ ->■ M + 7, 



(236) 



[37]. 
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pi -> mu + nu (1eV) at 5m 
pi -> mu + nu (0.2eV) at 5m 
K - mu + nu (1 .OeV) at 5m 
K -> mu + nu (0.2eV) at 5m 
-> e + gamma (1eV) at 100m 




size of wavepacket 

Fig. 9 The magnitudes of the probabilities in radiative and weak decays at rest measured 
at T = 0.33 x l(T 6 (/z), and T = 1.7 x 10- 8 (tt,K) [s] are shown. Parents and unobserved 
daughters are plane waves and observed particles are wave packets. The horizontal line 
shows the size of wave packets in units of crm 2 and the vertical line shows the enhancement 
factors at the finite distance, i.e., the ratios of the sum of the normal and diffraction compo- 
nents over the normal component. The decays of a pion to a muon and neutrino (red cross 
and green cross), and of a kaon to a muon and neutrino (blue-double cross and magenta 
box), and of a muon to a electron and photon (light-blue box) are presented. The masses 
are m v = 0.2 and 0.5 [eV/c 2 ], and rrCf f = 1.0 [eV/c 2 ]. 



in wide range. They have short mean life times and pseudo-Doppler effects may appear in 

M -> 2 7. (237) 

6.2. Finite-size correction 

The finite-size correction becomes large in the situation where the wave functions of the 
initial and final states overlap in wide area. This is realized at T < r and is important in 
slow decays of particles, such as weak decays and some gamma decays. Fig. (9) shows the 
enhancement factors at the finite distance, i.e., ratios of the total probabilities over the 
normal probabilities of the asymptotic region in various weak and radiative decays. For the 
large wave packets, the values become large. In this figure, initial states are plane waves and 
size of wave packet for the neutrino is expressed in the units of 1/tti 2 and is written in the 
horizontal line. The ratios p "™ a ' + ^ ,Jf are written in the vertical line. Pdiff 1S large in the 
region crm 2 > 10. Thus the finite-size corrections are non-negligible and important. 

In the region T 3> r, the finite size corrections vanish, and the decay rates are expressed 
by the standard formula. In this region, the number of the parents decrease as Noe~ T ^ T and 
that of the daughter becomes constant. 

6.2.1. Slow gamma decays of Nucleus . Photons produced from radioactive nucleus are 
measured through their interactions with nucleus in targets which have finite sizes. Hence 
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Fig. 10 The energy spectrum of photon in muon decays is shown at T = 1.7 x 10~ 8 [s] 
and m!y = 1.0 [eV/c 2 ]. The horizontal line shows the energy of photon in MeV at a^m^ = 
100, = o~ e = co and the vertical line shows the probability. The normal component (green) 
has a sharp peak around 54 [MeV], and the diffraction component(red) spreads over wide 
region in the lower energy side and resembles the background. 

S[T] expressed by wave packets describe the amplitudes of the process, 

N -»• N' + 7. (238) 

From Appendix C, the magnitude of light-cone singularity and the diffraction component 
are almost equivalent to those of point particles, and the total probabilities are modified by 

p diff- 

6.2.2. Muon decay to an electron and gamma. A muon decays to an electron and a photon, 

fi -»• e + 7, (239) 

where the photon's energy is about 50 [MeV], if the lepton number is violated. The lepton 
number violation has been observed in the neutrino oscillation phenomena but has not been 
in charged leptons. Precision measurements have been made and a new experiment has 
started [40]. Since the rate of this transition process is extremely small, it is important to 
know the corrections due to the pseudo-Doppler effect and finite-size correction. Those for 
the muon at rest of plane wave is studied here. From Fig. (6), the average energy of final 
state is larger than the initial energy. 

Fig. (9) reveals the enhancement of the rates due to the diffraction for plane waves of 
muon and electron and the wave packet for gamma. Fig. (10) shows the energy spectrum 
of 7 in the normal and diffraction components for <7 7 m 2 = 100, o" M = a c = 00 in the decay 
of muon of pn = at T = 1.7 x 10~ 8 [s]. The normal component has a sharp peak around 
E~i ~ 54 [MeV], whereas the diffraction component spreads over a wide region. Moreover 
the latter is much larger than the former in these parameters. Thus the corrections become 
important if the initial muon is a plane wave. The wave packet size of gamma can be 



"spec-diff-life-rest-scalar-T_5.0_m3_1 .0_m1_1 06_s3_1 0" + 
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Fig. 11 The energy spectrum of neutrino in pion decay at rest at T = 1.7 x 10~ 8 [s] is 
shown. The horizontal line shows the energy in MeV at a v m^ = 14.6, which corresponds 56 Fe 
and the vertical line shows the probability. Wave packets of another daughter and parent are 
oo and 0" paren tm^ = 10000. The neutrino mass is m u = 1.0 [eV/c 2 ]. The normal component 
(green) has a broad peak, and the diffraction component (red) spreads over in the low energy 
region. 

determined from the spectrum at the higher energy region of a known process, and is used 
for the calculation of the diffraction component of the present process. 

6. 2. 3. Weak decays . Neutrino measured through its interaction with a nucleus has a wave 
packet size of the nucleus. Hence the process of a nucleus 



is described by S[T]. The pion decay was discussed in the previous papers [41], and the 
neutrino's energy distribution is given in Fig. (11). From Appendix C, the magnitude of 
light-cone singularity and the diffraction component are about the same as point particles. 
The spectrum of diffraction component that gives the finite-size correction is distributed in 
the low energy region and that of the normal component is wide and has a peak. The peak 
is slightly shifted from that of the plane waves due to the pseudo-Doppler effect. From the 
shift and width of the normal component, the wave packet size can be determined and is 
used for the theoretical calculation of the diffraction component. 

6.3. Proton decay 

Proton is unstable and decays in Grand Unified Theory. In SU(5) GUT, a main decay mode 
is 




(240) 




(241) 
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The initial proton is in matter in ground experiments and final states are detected through 
wave packets. For large wave functions of a proton, neutral pion, and positron, they overlap 
in a wide area. For small wave functions, they overlap in small area. General cases with the 
symmetric wave packets 

a p , a n o, a e +, (242) 
of the four-dimensional momenta at positions 

(Pp,p p ;X p ,T p ), (p^ } p w o;X^o, TV), (p e+ , p e + ; X e+ , T e + ) , (243) 
are studied in the following. They are governed by an interaction Lagrangian 

%int = gipp(x)e(x)(p 7T (x). (244) 
The transition amplitude is an integral over (t, x) 

M (p^^ + e + ) =9 J dtj dxe-^ {3 - So)2 -^ {t - to)2 e R+ ^M, 

=g(2TTa s )l (2TTa t )h R+i *M, (245) 
for finite values of a s and at- M includes the spinors. a s and at are given in the expressions 

(246) 

Vr> V n a V P + \ /„._s 

— + — + —) . (247) 
a v a w o a e +J 

to and xo(t) are given in the form of Eq. (48) of an average velocity vo, 

f Up V w o V e+ \ 

Vq = cr s [ — -\ 1 . (248) 
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R and <f> in the exponent are obtained from Eqs.(50), (51), and (53) as, 



R — -^trajectory ~i~ ^momentum! (249) 
trajectory = " £ + 2<T S [ ]T ^ + 2<7 t ( £ 



where 



^momentum = \Ep{pp) ~ E n o(p n o) - E e +(p e +)j - ^-{p p ~ P> - p e +Y 



Pp=Pp-— (Pp ~ Pit ~Pe+), (250) 

p n o =p w o + — (p p -p w o -p e +), (251) 
<J e + 

Pe+ = Pe+ + — (Pp ~ Pn° ~Pe+), (252) 



and is a function of the momenta pj and positions Xj. 
1 Proton at rest 
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The proton in solid is at rest and is expressed with a small wave packet. In the rest system 
of the proton, the reduced momenta of Eq. (250) are 

Pp = — (Ptt° +Pe+), (253) 
p n o = p n o - —(p n o +p e +), (254) 

a e + 

p e + = p e + - — (p n » +p e +)- (255) 

a e + 

If the wave packet size of the positron is much larger than others 

a e + 3> a p ,a n o, v e + ~ v^o (256) 



then, 



a s = (257) 

= f + — 7 X ;{Ve+ ■ VTfl) 

a t a p + a n o a e + a e +{a p + a 7T o) 
„2 

(258) 



dp + <Jtt 



we have 

^momentum = ~y (^p(fp) ~ E n a(p n a) - E e +(p e+ )^j - y (j?^ + Pe+f ■> (259) 

Pp = —{p n o +Pe+), Ptt° =Ptt«, Pe+ = Pe+ ■ 
Op 

Thus in a region of large a s and at, the conservation law of the momentum and energy are 
the same as that of plane waves, but in small a s the momenta are spread over in wide region 
and the energy conservation law is modified. For the large at, in the event of 

Piro+Pe+^O, (260) 

the energy conservation takes the form, 

Ep(pp) ~ £*o(^ro) - E e+ fn 0. (261) 

p p could be very different from p p = 0, hence the modified conservation law derived from 
pseudo-Doppler effect should be taken into account for the experimental analysis in this 
region. 

Since at is finite, the decay probability is proportional to T in the region T <C T proton , 
and the decay rate is constant in a wide range of T despite the fact that the spectrum 
is distorted, where T pro ton is the average life-time. Thus the proton at rest decay with the 
constant rate even at small T and the proton decay experiment is doable if the life-time is 
less than 10 34 — 10 35 years. 

6.3.1. Other decay processes . Three body decays such as /i — > e + v + n — > p + e + v 
and others have light particles in the final states and are modified by the pseudo-Doppler 
and finite-size corrections. They will be presented in a separate paper. 
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6.4- Thermodynamics of small quantum particles 

When an excited state of heavy atom of the large wage packet size makes transition without 
changing the momentum and emit a photon, they follow the modified energy conservation 
law. If the atoms are in thermodynamic equilibrium of a temperature T, the state of the 
energy E follows a distribution 

p(E,p)=N(-PE), (262) 

where j3 is inversely proportional to the temperature, and N(—f3E) becomes Planck 
distribution for Bosons and Fermi-Dirac distribution for Fermion. 

In the situation where the wave packet size of the atoms is much larger than the wave 
packet size of a photon and the atoms are bound together strongly similar to Mossbaur 
effect, a distribution of the photon receives pseudo-Doppler effect and Mossbaur like effect. 
Then the temperature of photons becomes to deviate from that of the atoms. 

From Eq. (57), in a situation 

<7 7 <C a a, (263) 

the velocity vq agrees with the velocity of the photon. The photon's energy is given by 
Eq. (208), hence the energy distribution of the photon emitted from atoms is 

p(R y ,P) = N-J—-J = ^. (264) 
e PE~, _ i k 

Thus the effective temperature of the photon is k times of that of the atoms. 
7. Summary and implications. 

We have developed a theory of the diffraction induced by many-body interactions and com- 
puted the finite-size corrections to the rates of slow transitions caused by electromagnetic 
and weak interactions. The large corrections to the Fermi's golden rule Eq. (13) were found 
in certain processes. 

The Fermi's golden rule is applicable to the rates in the particle zone where the initial 
and final states separate completely and their wave functions do not overlap. The rates are 
not subject to 1/T correction. In the wave zone, however, the 1/T corrections are found 
using the wave functions that satisfy the boundary conditions. Because they have universal 
properties, they are observable in scattering experiments. The finite-size correction reveals 
the diffraction pattern of one quantum interference. The intermediate-time region of particle 
decays when the parent and daughters co-exit with a finite overlap of wave functions is the 
example of wave zone. The state is a superposition of the parent and daughters and has a 
finite expectation value of Hi nt . Thus the kinetic energy varies here, the decay rates and 
scattering cross sections are different from their asymptotic values. 

The finite-size corrections are inevitable consequences of the boundary conditions at T 
and have the magnitude that depends on the sizes of a s . The size of a s can be controlled, 
and the finite-size corrections will be verified in experiments. Especially if the coherence 
length, hE/(m 2 c 3 ), is a macroscopic size much larger than de Broglie wave length h/\p\, 
they would reveal in macroscopic scales. For neutrinos or photon of the effective mass of the 
order [eV/c 2 ], hE/(m 2 c 3 ) becomes a macroscopic size. Hence the finite-size correction may 
become visible in a macroscopic distance. 
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If a s = finite, at = finite, the wave functions overlap only in microscopic regions. Waves in 
experiments at macroscopic distance are in the particle zone, and the finite-size correction 
disappears. Nevertheless, the probability receives pseudo-Doppler effect even in this region. 
The distortion of the energy distribution becomes stronger with smaller wave functions and 
may become visible. It becomes drastic if both of Mossbaur and pseudo-Doppler effects are 
combined as in Eq. (208). The final state of huge kinetic energy, much larger than the initial 
one, is formed with a small probability. The conservation of kinetic energy is violated in each 
event but satisfied for the average value over the classical time interval. Despite the fact that 
S[oo] conserves the kinetic energy, the modified energy 5E conserves approximately. In fact, 
these behaviors would have been considered as artifacts of the detectors and absorbed in cal- 
ibrations of the detector. The present results might help toward the complete understanding 
in this direction. 

In slow radiative and weak decays of particles A! and N' of plane waves, 

A' -»• A + 7, N' -> N + v, (265) 
the rates and energy distributions of 7 or v in the asymptotic region, 

r total = r(°>, (266) 

Vtotai(E) = V^(E), (267) 

are computed with plane waves and ie prescription, where E is the energy of the observed 
particle. Our results for 7 or v measured at L reveal the finite-size corrections and pseudo- 
Doppler effects and are expressed in the form, 

T total = rW+rWfl(L;<r), (268) 
Vtotai(E) = V^(E;a)+V (dlff \L,E;a). (269) 

The diffraction components have magnitudes summarized in Eqs. (154) and (178), and are 
important in various processes of leptons, hadrons, nucleus, and positoronium and negligible 
in ordinary atoms. At L 3> Lq 

rWfl(L;(j)4 0, (270) 
p( di ff)(L,E;a)^0, (271) 
where Lq is the minimum value of the mean length and coherence length, 



HE 

min <j cr, — ^-3 \ ■ (272) 



As a -)• 0, 



\V {n) {E-a) -V {G) {E)\ -»• large, (273) 

T^ dlff) (L; a), V (diff) (L, E; a) -> 0, (274) 

whereas as a — > large, 

\V {n) {E;a) -V {0) (E)\ ^ 0, (275) 

(L; cr), V {diff) (L, E; cr) — >■ large. (276) 

The normal and diffraction components behave differently with cr, and Figs. (6) and (9) show 
them. From Eqs. (273) and (275), a gives the effects to observables regardless of its magnitude 
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and can be determined experimentally from the energy spectrum Figs. (5)-(8) of the normal 
components V^ n '(E). The values computed by Fermi's golden rule have large corrections, 
Figs. (5)-(ll), and the theoretical values can not be consistent with the experiments without 
corrections. Hence it would be easy to verify the finite-size corrections or pseudo-Doppler 
effect. Because the Fermi's golden rule is applied to many problems in wide area, it would 
be important to confirm the corrections. 

The decay rate of proton in matter is constant in T< T pro t on , if the final states are 
measured due to the boundary conditions. This agrees with the standard one, and the 
proton decay will be detected if GUT is correct. Finally unusual luminescence and thermo- 
dynamic properties of quantum particles caused by overlap of wave functions will be verified 
in experiments. 

Constituent particles such as molecules, atoms, nucleus, or elementary particles have small 
intrinsic sizes and are expressed with wave functions of finite sizes in certain situations. 
Consequently, their reactions may be affected with finite-size corrections or pseudo-Doppler 
effects, even though any measurement are not made. Physical system may show such unusual 
behavior as the non-conservation of kinetic energy. Nevertheless the average energy over 
a long period recovers the conservation law. Hence the phenomena may appear in non- 
stationary and time-dependent processes. [33-35]. Macroscopic quantum phenomena in this 
situation have been barely studied and will be discussed in subsequent works. 
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Appendix A. Finite-size correction to Fermi's golden rule 

A -I. Approximation with Dime's delta function 
Integrals over the finite interval 

r d te^ = e^ S ' m{UjT/2 \ (Al) 
Jo w /2 

j\ ldh e^~^ = (^^f , (A2) 
are approximated normally with 



/ dte iujt = 2tt5(uj), (A3) 
Jo 

[ dhdt2e iuj{tl - t2) = 2ttT5(oj), (A4) 
Jo 



for a large T. They have been applied in computing a decay rate and cross section and are 
explained in most textbooks. 
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Finite-size correction to this formula depends on oj. If co is discrete, 

T; u = 0, 



dte 



iuit 



2sin(wT/2) JujT/2 _ 



dt 1 dt 2 e i ^ tl - t2 



T ; u = 0, 



2sin(wT/2) 



; u^Q, 



and the averages over finite time interval 5T of 5Tuj S> 1 are 



Aver 



Aver 



dte 



iwt 



T; w = 0, 



dtxd^e^-^ 



T , W = 0, 



For the average probability at a finite-T, the correction is given by 

2 



CO 



2^2' 



(A5) 



(A6) 



(A7) 
(A8) 

(A9) 



A-II. Correction by Taylor expansion 

If co is continuous, there exists states of infinitesimal energy differences. The correction 
becomes non-trivial and is studied here. The following integral for co\ < < co 2 



I(ui,u} 2 ;T) = / dcog(co) 



sin(wT/2) 



co 



(A10) 



coincides with 2irTg(0) if the second equation of Eq. (A3) is used. To find next order terms 
in 1/T, we expand g(co) 



g(u) = g(0) + £ 
and change the variable to x = wT, 



i=i 



g (0 (Q ), 
il 



-co 



(All) 



I(coi,co 2 ;T) =T / dx 

Ju> 1 T 



sin(x/2) 



g(x/T) 



v-^ f (0 (0) r 2T feun(x/2y 



dx 



(A12) 

l 

The integrand in the above equation is finite at x = for 1 = 0, but those for I > 1 vanish. 
At large x, sin 2 (x/2) « 1/2. So we have 



/K, W2 ;T)«2^(0) + ^^T f 

Z!T J <* 



i 1 1 

x 2 2 



2ttT 5 (0) + ^— ^log 



COo 



Z>2 



2/! 



i-1 



(A13) 
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Choosing uj\ = —uj 2 , we have 

I(-u, 2 ,u, 2 ;T) = 2vrT 5 (0) + £ f^^j(^l(l - (")')), (A14) 



l> 



The second term in the above equation is the 1/T correction. This correction depends on 
the cut off frequencies and the constant term i.e., uj 2 term vanishes at the symmetric cut off 

oj\ = —co 2 - So lim I(—uj 2 ,uj 2 ;T) agrees with 2irTg(0). The finite-size correction is written 

cj 2 — >o 

in the form, 

I(-oj 2 , cj 2 ; T) = 2vrT<7(0) ( X + ^) ' (A15) 
where To is roughly estimated with a size of atom as 

io~ 10 M 1Hri fk , 

T = — = 0.3 x 1(T 18 s , A16 

c 

in atomic physics. Hence in an experiment of the size 1 [m], T = 1/c = 0.3 x 10~ 8 [s], and 
the correction becomes 

Tn 



= io- 10 

r 

This value is negligible and the finite-size correction vanishes at a macroscopic distance. 



= 10~ iU . (A17) 



Appendix B. Level density and correlation function :quantum mechanics 

We summarize the probability at a finite T, Eq. (11), of systems of various level densities. 
When a level density p(E) is given, a number of states below E, s, satisfies 

The correlation function is expressed with s in the form, 

f'OO f'OO 

g (+)( tl -t 2 )= dEp(E)e^ E - E "^^ = ds j(E(s)-E )( tl -t 2 )^ (B2) 
Je Jo 

g (-)( tl - t 2 ) = [ E ° dEp(E)e^ E - E "^-^ = f° dse< E ^- E ^-^ . (B3) 

JE m Js m 
Using s, we have the integral over a finite interval of s of Eq. (A3), C(T) = 1(0, oo,T), 

C(T) = I dt l dt 2 g(t 1 -t 2 ) = I ' 2 ds I ^i^ 2 e Ms)(tl "* 2) , (B4) 

J0 J s t Jo 

w(s) = E(s) - E , g(h - t 2 ) = g(+\t! - t 2 ) + g^(h - t 2 ), 

which agrees with 

= 2ttT / ds5(u)(s)) = 2vrT— — — = 2irTp(E ), (B5) 
J Sl K(*o)| 

if oj(s) = has a simple root, so in s± < sq < s 2 and uj'(sq) is not too small. This is equivalent 
to Eq. (^43). C(T) is written also in the expression 

C(T) = T / T dfr(f) + f° d^g(0 - C (B6) 

J-T J-T JO 
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B-I. Regular level density 

For the level densities that are regular at E = Eq, 



p(E) = c , 



1 



E- Eo-iT 1 



-c\E-E \ 2 



(B7) 



we have 



Table Bl 



p(E) 


S 


g(h - t 2 ) 


C(T) 


co 


c {E - E ) 


2irc 5(ti - t 2 ) 


27TCqT 


1 


J 


r* dE' 


e -T(t t -t a ) 




(' 






E - Eo-iT 


j En E'-Eo- iT 


TT j 


e -c\E-E \ 


J 


( ' E dE'e- c \ E '- E °\ 2 

Eq 


(il-*2) 2 

e 2c 


\ 


Y 







The correlation functions g(t\ — t 2 ) in Table Bl are short range and Eq. (A3) is applicable. 
C(T) are proportional to T and the corrections are proportional to 1/T. 

B-II. Weakly singular level density 
For the level densities, 



Table B2 



P{E) 


P(0) 


P'(0) 


s 


g{h - 1 2 ) 


C(T) 


e -C\E-E \ 


1 


± 


C dE'e- c \ E '- E °\ 
Je„ 


l 

c - i{h - 1 2 ) 


^T7r-2Clog^ 


C(E - EoY 





CO 


C{E - E ) p+1 


r(£±i)2 p -i 


Txoo 


(0 <p< 1) 


p+l 


r(-f)|*i-tari 



The correlation functions g(ti — £2) are long range and Eq. (.A3) is applicable in the first one 
but is not in the second one. C(T) are proportional to T and the correction is proportional 
to (logT)/T in the former one and the proportional constant diverges in the latter one. 
Especially in the latter case, the level density satisfies, p(Eq) = 0, but the coefficient C(T)/T 
diverges. 

B-III. Singular level density 
For the singular level densities 

Thus Eq. (A3) is not applicable in the level densities of Table B3. The level density satisfies 
p(Eq) = 00 and decay rates C(T)/T is proportional to a positive power of T. Thus this 
system should reveals unusual non-Markov property. 



54/57 



Table B3 



p{E) 


P(0) 


P'(0) 


s 


g{h - h) 


C(T) 


C5{E - E ) 


oo 


oo 


C9(E - Eo) 


c 


CT 2 


— °^ , ; <p < 1 


oo 


oo 


c(e - E y- p 


[2pwT(l-p) 


2T i+ P 


\E-Eq\p 


l-p 


"V 7T 2|*i -t 2 \ l ~P 


P 



B-IV. Relativistic wave packet 

The spectrum density for a two-body decay ^4 — > B + C of relativistic particles Eq. (128) 
integrated over pc is 

PH = / ^e-^-v°-^ 2 5(u - E A (p A ) + E c (po) + E B (p A - p c ))- (B8) 

p(0) is finite because the root of ui = exists and /o(w) satisfies 

JO, cj>^a, mQ . 
[e -i , u) — >■ — oo. 

Due to smooth asymptotic behavior Eq. (-B9), the following integral converges even at a 
finite T, 

/ <b (^ffl)'^) = T (2*0) / d. (%2) 2 + i,) , (BIO) 



where 77 is given by 



' du ( Sm(tjT/2) ^ p(u) > 0, = p(u>) - p(0), (Bll) 



Co' 



and the positivity Eq. (35) is fulfilled. 



Appendix C. Light-cone singularity for general composite systems 

C-I. Light-cone singularity and form factor 

Composite particles such as hadrons, nucleus, atoms, and molecules have internal structures 
and have modified light-cone singularities. The form factor F((p A — Pc) 2 ) m 

(C;p c \J(x)\A;p A ) = e-^-^- x rF((p A -p c ) 2 ), (CI) 

depends on Lorenz scalar (p A — pc) 2 , where J(x) is source operator of 7, v, or others that 
are detected. Lorentz indices are ignored. The correlation function is 

AA ' c(5x) = ji^ J ^^\ F ^-Pcf)\ 2G (PA,Pc)e- i{pA - pc) - Sx 

G(p A ,p c ) = Y,( rr *)A,c, (C2) 



55/57 



and we have the light-cone singularity, 



(2nf 



— |i^-pc) 2 )| 2 e-^-^)- fa 



E(pc) 
2 

i7r(2-7r) J 



d qlm 



_q 2 + m 2 A - m A c + 2pA ■ 



q — le 



2^ 1 2 ^—iq-5x 



\F{m z c - m\)\ z —8{X)e{5t) + "others" 



(C3) 



Thus we have 



A A , c (5x) = G(p A , -iJL- PA )\F(m 2 c - m 2 A )\ 2 (±-5{\)e{5t) + "others" 



The magnitude is renormalized by the form factor \F(m%-, — m A )\ and the form is kept 
intact. 

C-II. Strength 

\F(mQ — rn A )\ 2 of various systems, using the energy gap, AE, and size, R, 



m A =m c + SE, F(q 2 ) = F(0)e~- 



(C4) 



is written as 



F(m 2 c - m\) 2 = F(0) 2 e- R ^^, 



(C5) 



and is determined by R\/2mcSE. 

The typical values for bound states composed of electron, nucleon, quark, and //-on, eA, 
e + e~ , NN, qq, qqq, fiN, fi + e~, and e(K)A are 



R = < 



— 2~ ; (atom: eA), 



rUeC* a 
he 1 



m P .cr a 

h 



— ; (positronium : e e 



h 

m„c 



; (nucleus : NN), 
; (hadron : qq, qqq), 



he l , . 

tt — ; {/j, atom: fiA), 

n^c 1 a 

he 1 . . 
— 2~ ; \n atom: fie), 



m e cr a 
he 1 
, m e c 2 Nkoi 



; (K-shell atom: eA), 



AE = < 



m e c 2 a 2 



; (atom: eA), 



2^,2 



m e c a 
2 

m^c 2 
100 

V 



■; (positronium : e + e 



(nucleus : NN), 



m q c 2 ; (hadron : qq, qqq), 



2„,2 



m P c a 



; (// atom: fiA), 
■; {fjL atom: /xe), 



2^,2 



m P c a 



■; (K-shell atom: eA), 



(C6) 
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and, 



m c 



rriN] (atom: eA), 

2m e ; (positronium : e + e~), 

Attin; (nucleus : NN), 

Bm q c 2 ; (hadron : qq, qqq) , 

?7i TVj (a* atom: fJ,A), 

m M ; (fj, atom: fie), 

mjsr; (K-shell atom: eA). 



(C7) 



We have 



i? 



M 2 c 2 




mAr- 



m e c 2 a 2 



m^ 



M 2 2 

1; (positronium : e + e~ 
< 1; (light nucleus), 

lOOm^ 

w3(4= 100); (nucleus : 2ViV), 
« 1; (hadron : qq, qqq), 
3; (/i atom: fiA), 



> 50; (atom: eA), 



(C8) 



14; (/x atom: /ie), 

- = 10; (K-shell atom: e^4), 



and 



2 \2 



F{mh - ml) 



F(0) 2 x 0(1); hadron, positoronium, light nucleus, 
F(0) 2 x O(10 _1 ); /iiV-atom, heavy nucleus, 
F(0) 2 x O(10~ 5 ); ixe-atom, K-electron, 
F(0) 2 x O(10- 10 ); atom. 



(C9) 



Thus the magnitude of light-cone singularity is about 1 of F(0) 2 for positoronium, light 
nucleus, and hadrons, and 10 _1 for /iiV-atom and heavy nucleus, and 10 -5 for tie-atom and 
K-electron, and 10 -10 for atom. Atom is composed of heavy nucleus and electrons and its 
size and energy gap are determined by electron's mass. Hence the large ratio ^/m^ /m e 
determines the overlap of an atom C with an atom A, and the strength of the light-cone 
singularity becomes extremely weak, accordingly. For other particles composed of equal 



masses, F(m\ 



F(0). 
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